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Abstract 


Aircraft  collisions  with  avian  species  are  a  serious  safety  problem  as  well  as  a  serious 
economic  issue.  Aircraft  /  bird  strikes  have  resulted  in  33  fatalities,  the  loss  of  39  aircraft,  and 
damages  to  aircraft  in  excess  of  $820M  for  the  United  States  Air  Force.  The  objective  of  this 
paper  is  to  create  a  closed  form  mathematical  model  that  estimates  the  probability  of  a  bird  / 
aircraft  collision  and  provides  a  risk  score  that  can  be  utilized  to  underpin  decisions  made  by 
planners  and  pilots. 

The  major  components  of  the  model  are  the  spatial  Poisson  process,  the  extended  spatial 
Poisson  process,  a  gamma  distribution  of  bird  altitudes,  a  relative  risk  score,  a  standardized  risk 
score  scale,  and  a  risk  filtering  and  ranking  method.  The  spatial  Poisson  process  allows  for  an 
independent  distribution  of  birds  within  a  bounded  area.  The  extended  spatial  Poisson  process 
accounts  for  the  removal  of  birds  from  calculations  within  the  bounded  area  after  they  have  been 
encountered.  The  gamma  distribution  models  the  distribution  of  specific  bird  altitude  bands 
within  a  bounded  area.  The  relative  risk  score  is  a  weighted  risk  score  for  19  different  species  of 
birds  that  an  aircraft  might  encounter.  The  standardized  scale  aggregates  all  risk  scores  over  all 
the  bird  species  and  then  calculates  the  value  in  a  0  to  10  scale.  The  risk  filtering  and  ranking 
model  combines  the  effects  of  a  hit  with  the  likelihood  of  a  hit  and  displays  the  result  in  a 
graphic.  The  overall  model  that  combines  these  components  and  calculates  the  output  is  an 
original  contribution  to  the  field  of  aircraft  /  avian  collision  models. 

Exercising  the  model  reveals  significant  factors  that  influence  the  risk  score  associated 
with  flying  in  a  particular  area.  They  are  the  total  number  of  birds  in  the  bounded  region,  the 
mix  of  species  within  the  bounded  region,  the  size  of  the  aircraft,  and  the  gamma  height 
distribution  of  the  birds  within  the  bounded  region. 
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Knowing  the  gamma  height  distribution  for  the  specific  birds  in  an  operations  area  (AO) 
can  provide  more  fidelity  to  the  planner.  In  fact,  in  several  scenarios  where  the  same  number 
and  species  of  birds  for  an  AO  was  used,  the  difference  in  the  overall  aggregated  risk  score  was 
twice  as  high  as  the  score  that  was  calculated  when  the  gamma  height  distribution  was  not 
known.  Additionally,  when  there  were  densely  populated  altitude  bands  of  birds  in  the 
operations  area,  avoiding  these  bands  cut  the  overall  risk  score  by  up  to  50%.  This  is  very  useful 
information  for  decision  makers  to  have  when  they  are  planning  the  specifics  of  their  operations. 
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ESTIMATING  BIRD  /  AIRCRAFT  COLLISION  PROBABILITIES  AND  RISK 
UTILIZING  SPATIAL  POISSON  PROCESSES 

1.  Introduction 

1.1.  Background 

Random  collisions  of  relatively  small  numbers  of  entities  in  very  large,  yet  bounded, 
spaces  are  rare,  but  not  impossible.  In  fact,  consider  the  time  period  spanning  15  January  2009 
to  4  February  2009.  During  this  20  day  period,  the  U.S.  commercial  Iridium  spacecraft  was  hit 
by  a  defunct  and  out  of  control  Russian  satellite  at  an  altitude  of  about  800km  (500  miles)  over 
Siberia  (BBC  NEWS  2009);  additionally,  a  Royal  Navy  nuclear  submarine  was  involved  in  a 
collision  with  a  French  nuclear  sub  in  the  middle  of  the  Atlantic  Ocean  (BBC  NEWS  2009)  and 
finally,  US  Airways  flight  1549  made  an  emergency  water  landing  as  a  result  of  a  bird  strike  (a 
large  flock  of  Canada  geese  disabled  both  engines)  in  the  Hudson  River  in  New  York  (U.S. 
Airways  2009). 

While  the  first  two  cases  of  collisions  are  quite  an  anomaly,  aircraft/wildlife  strikes  in 
general  do  happen  frequently.  According  to  the  Federal  Aviation  Administration  (FAA),  which 
has  maintained  a  wildlife  strike  database  since  1990,  there  have  been  121,000  (civil  and  U.S.  Air 
Force)  wildlife  strikes  between  1990  and  2010.  The  U.S.  Air  Force  (USAF)  has  been  collecting 
statistics  since  1985  and  has  recorded  95,383  wildlife  strikes  spanning  the  period  from  1985  to 
2011,  resulting  in:  33  fatalities,  the  loss  of  39  aircraft,  and  damages  to  aircraft  in  excess  of 
$820M  (U.S.  Air  Force  Safety  Center  2012). 

According  to  the  FAA,  92%  of  the  bird  strikes  to  commercial  aircraft  occur  at  or  below 
3,500  ft  above  ground  level  (AGL).  Using  this  same  threshold,  the  USAF  has  recorded  96.72% 
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of  its  bird  strikes  below  3,500  ft  AGL.  In  order  to  narrow  down  scope  of  the  problem  set, 
statistics  on  USAF  wildlife  strikes  by  phase  of  operations  were  examined.  The  following  table 
highlights  the  number  of  bird  strikes  for  the  different  phases  of  operation  (U.S.  Air  Force  Safety 
Center  2012). 


Table  1  USAF  Wildlife  Strikes  by  Phase  of  Operations  (1995-2011) 


Phase  of  FI ight 

Cost  ($MM) 

%  of  Total 

Count 

%  of  Tota 1 

Takeoff/Initial  Climb 

137.0 

32.2% 

7299 

12.3% 

Enroute/Ai  r  Work/Ai  r-to-Ai  r/Ai  r  Refue  1  i  ng 

19.2 

4.5% 

2830 

4.8% 

Flight  Demonstration 

1.9 

0.4% 

27 

0.1% 

Low  Level/Air-to-Ground/Air  Delivery 

174.2 

40.9% 

6949 

11.7% 

FI  over 

0.0 

0.0% 

10 

0.0% 

Traffic  Pattern/Go-Around 

31.4 

7.4% 

7618 

12.8% 

Initial  Approach/Final  Approach/Landing 

46.8 

11.0% 

16048 

26.9% 

Parked/Ground  Ops 

3.6 

0.8% 

423 

0.7% 

Unknown 

12.1 

2.8% 

18400 

30.9% 

These  statistics  clearly  do  not  support  the  “big  sky,  little  plane”  theory  and  as  such,  there 
are  extensive  mitigation  and  alert/advisory  models  and  systems  in  place  throughout  the  world  to 
help  reduce  the  occurrence  of  aircraft/bird  collisions.  Both  active  and  passive  mitigation  & 
abatement  methods/systems  are  employed  around  airfields  in  order  to  diminish  the  possibility  of 
an  aircraft/bird  strike.  These  methods  include:  pyrotechnics,  bioacoustics,  exclusion,  propane 
gas  cannons,  falconry,  dogs,  radio  controlled  crafts,  grass  habitat  management,  sanitary  landfill 
management,  agricultural  leases,  and  depredation  (U.S.  Air  Force  2004). 

However,  for  most  low  level  routes,  air-to-ground,  and  air  delivery  operations,  mitigation 
and  abatement  methods  and/or  systems  are  not  available  and  as  Table  1  indicates,  the  largest 
dollar  cost  and  a  high  total  count  is  evident,  even  though  a  smaller  percentage  of  flights  fly  this 
profile.  To  assist  in  reducing  the  probability  of  a  bird  strike  in  this  phase  of  operation,  advisory 
models  and  systems  are  in  place  to  assist  with  planning  and  real-time  avoidance.  These  advisory 


2 


and  alert  systems  vary  on  methods  used  (historical,  real-time,  algorithms,  etc)  and  provide  pilots 
and  planners  with  different  levels  of  infonnation  and  fidelity. 

1.2.  Problem  Statement 

The  purpose  of  this  research  is  to  create  a  closed  form  mathematical  model  to  detennine 
the  probability  of  an  aircraft/wildlife  collision  encounter  (strike)  in  a  defined  space  (e.g.  on  a 
segment  of  a  Low  Level  route/Air-to-Ground/Air  Delivery)  with  a  set  of  given  parameters  and 
then  give  weighted  risk  associated  with  the  encounter.  The  output  analysis  will  be  presented  via 
a  Decision  Support  System  (DSS)  in  order  to  underpin  decisions  by  pilots  and  planners  when 
coupled  with  existing  alert  and  advisory  systems. 

1.3.  Research  Objectives 

To  detennine  the  probability  of  an  aircraft  bird  encounter  that  has  the  possibility  of  leading  to 
a  bird  strike,  the  following  research  objectives  are  considered: 

•  Study  the  systems  that  provide  avian  detection,  alerts  and  advisories;  discuss  the  level  of 
fidelity  they  provide. 

•  Detennine  avian  characteristics,  geographical  and  environmental  conditions,  and  other 
factors  that  drive  avian  population  densities  and  spatial  patterns. 

•  Conduct  research  on  Spatial  Poisson  Processes  and  relate  them  to  aircraft/avian 
encounters  for  a  defined  space. 

•  Highlight  and  discuss  existing  aircraft/avian  collision  mathematical  models. 

•  Conduct  research  on  Extended  Spatial  Poisson  Processes  and  how  they  relate  to  the 
encounter  model  and  techniques  for  solving  matrix  exponentiation. 
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2.  Literature  Review 


2.1.  Overview 

This  section  discusses  the  completed  research  that  support  the  research  objectives  listed 
above. 

There  is  a  vast  amount  of  research  that  has  been  completed  on  bird  strikes.  In  fact  there 
are  numerous  international  committees  that  meet  to  discuss  trends,  best  practices,  research, 
detection  and  advisory  models,  and  other  related  topics.  Additionally,  there  are  many  books, 
journals  and  newsletters  that  have  been  published  which  discuss  this  issue.  Finally,  there  has 
been  a  marked  increase  in  the  technology  developed  for  avian  detection  and  advisory.  All  of 
this  is  aimed  at  reducing  the  number  and  severity  of  bird  strikes  throughout  the  world  as  both  the 
number  of  aircraft  operations  and  large  bird  species  populations  increases.  For  example,  in 
North  America,  13  of  the  14  largest  (>3.6  kg  body  mass)  bird  species  have  shown  significant 
population  increases  in  the  past  40  years.  In  fact,  the  migratory  and  non-migratory  population  of 
Canadian  geese  (average  weight  of  4.2  kg)  has  more  than  quadrupled  from  1.2  million  to  5.5 
million  birds  in  North  America  from  1970  to  2008.  (Dolbeer  2009). 

Random  collisions  in  Md-  space  have  been  studied  at  great  length  by  many  scholars. 
There  are  numerous  areas  of  specific  study  within  the  three  main  groups  of  models  which  are: 
gas  particle  models,  satellite  models  and  historical  models.  This  paper  will  look  at  some  of  the 
specific  areas  of  study  related  to  random  collision.  The  focus  will  be  on  the  Spatial  Poisson 
Process,  as  it  allows  the  most  flexibility  when  dealing  with  the  characteristics  and  behaviors  of 
aircraft  and  bird  interactions. 

The  end  result  of  the  research  is  to  provide  a  Decision  Support  System  that  combines 
some  of  the  ideas  garnered  from  the  review  of  avian  characteristics,  advisory  &  alert  systems 
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with  the  mathematical  modeling  of  spatial  Poisson  processes.  The  goal  is  to  create  a  model  that 
allows  for  parameter  inputs  for  1)  the  area  of  operation  (AO),  2)  the  type,  number  (and  other 
associated  characteristics)  of  the  birds  in  the  AO  and  3)  the  size  (wingspan)  and  distance  of  the 
aircraft  flying  through  the  AO  and  then  provides  a  probability  for  the  aircraft  /  bird  encounter 
along  with  a  risk  value. 

2.2.  Avian  Characteristics  and  Environmental  Factors  that  Influence  Populations 

In  order  to  derive  avian  density  for  a  particular  area,  one  must  first  discuss  some 
fundamentals  of  avian  movement  patterns.  Avian  movements  can  be  partitioned  into  three 
categories:  migrating  birds,  commuting  birds  and  resident  birds.  (FAA  2010).  Within  these 
categories,  there  are  many  factors  that  affect  the  avian  population  density  for  a  specific  area.  The 
most  influential  factors  affecting  bird  movements  and  roosting,  loafing  and  feeding  locations  are 
wetlands  locations,  land  use  practices,  water  management  facilities  and  agricultural  activities 
(Cleary  and  Dolbeer  2005).  Specifically,  the  most  influential  areas  within  each  group  can  further 
be  defined. 

Land  use  practices  that  attract  birds  are: 

•  Waste  disposal  operations 

•  Underwater  waste  discharges 

•  Trash  Transfer  Stations 

•  Composting  operations 

•  Fly  ash  disposal 

•  Recycling  centers 

•  Constructions  and  demolition  debris  facilities 
Water  management  facilities  that  attract  birds  are: 
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•  Storm  water  management  facilities 

•  Wastewater  treatment  facilities 

•  Artificial  marshes 

•  Wastewater  discharge  and  sludge  disposal 
Agricultural  activities  that  attract  birds  are: 

•  Crop  production 

•  Livestock  production 

•  Aquaculture 

Therefore,  it  is  noted  that  the  population  density  of  birds  varies  considerably  among 
regions  and  habitats.  In  general,  greater  numbers  of  bird  species  are  attracted  to  areas  offering 
more  diverse  food  sources;  an  abundance  of  food,  cover  and  water  leads  to  larger  numbers  of 
birds. 

The  distribution  and  density  of  birds  also  changes  with  the  season.  In  the  Northern 
Hemisphere,  bird  numbers  peak  in  the  summer,  after  the  breeding  season.  In  the  far  north,  Snow 
Geese  breed  in  colonies  of  up  to  150,000  pairs.  Large  sea-bird  colonies  comprised  of  thousands 
of  nesting  birds  are  found  along  both  the  eastern  and  western  seaboards  of  Canada  and  the 
United  States.  Around  the  Great  Lakes,  on  small  islands,  gull  colonies  have  been  documented  as 
containing  over  40,000  breeding  pairs.  During  migration,  birds  of  some  species  funnel  to  and 
congregate  at  key  staging  areas  along  the  flyways.  As  a  result,  relatively  small  areas  can  become 
the  temporary  home  to  extremely  high  concentrations  of  birds,  and  many  airports,  low-level 
routes,  operating  areas  and  ranges  are  located  along  major  migratory  bird  routes  (Transport 
Canada  2004).  The  migratory  phenomenon  is  worldwide.  Large  eastern  European  birds  migrate 
to  Africa  annually,  passing  along  the  eastern  shore  of  the  Mediterranean  Sea.  Additional 
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European-Africa  routes  are  over  the  Iberian  Peninsula  and  the  Italian  Peninsula  to  central, 
southern  and  eastern  Africa  (Eschenfelder  2005). 

Additionally,  the  time  of  day  affects  the  population  activity.  The  highest  levels  of  daily 
wildlife  activity  normally  occur  +/-  one  hour  of  sunrise/sunset  as  birds  move  to  and  from  their 
roosts  (U.S.  Air  Force  2011). 

Finally,  the  most  significant  influence  on  the  altitude  at  which  birds  fly  is  the  weather, 
specifically  cloud  cover  and  wind  fields.  Birds  may  fly  lower  when  it  is  cloudy  or,  if  the 
overcast  is  not  too  thick,  they  may  ascend  through  it  to  reach  the  clear  skies  above.  If  favorable 
tail  winds  are  to  be  found  in  certain  altitudinal  strata,  birds  often  ascend  or  descend  in  order  to 
take  advantage  of  them.  A  comparative  analysis  of  the  influence  of  weather  on  the  flight 
altitudes  of  birds  was  conducted  by  researchers.  They  found  that  flight  altitudes  of  birds  differ 
among  species  and  vary  greatly  from  day  to  day.  Many  factors,  in  addition  to  weather,  may 
influence  the  flight  altitudes  of  birds.  The  study  utilized  several  different  parameter  inputs  and 
regression  analysis  was  accomplished  for  different  species  of  birds  in  order  to  form  a  prediction 
model  for  bird  altitude.  Data  was  collected  via  radar  to  the  level  of  individual  birds  and  type  was 
determined  by  visual  or  classification  of  wing  beat  frequency.  Some  of  the  input  parameters 
were:  relative  humidity,  min  and  max  daily  temperature,  lift  index,  sea  level  pressure,  wind 
speed,  aeronautic  index,  boundary  layer  height,  thermal  index,  vertical  wind  gradient  and  cloud 
cover.  Different  combinations  of  these  factors  explained  40%  to  70%  of  the  variance  in 
maximum  flight  altitudes  for  the  particular  species  studied  (Shamoun-Baranes,  Van  Loon,  et  al. 
2006). 


7 


2.3.  Avian  Detection  and  Advisory 


Different  techniques  have  been  used  to  warn  the  pilots  and  flight  schedulers/planners  of 
potential  bird  threats. 

Radar  has  been  used  extensively  to  monitor  bird  movements  and  warn  the  relevant 
personnel  (Sodhi  2002).  One  of  the  primary  radars  used  to  monitor  bird  movements  in  the  U.S. 
is  the  WSR-88D  NEXRAD  radar.  The  WSR-88R  is  the  Doppler  weather  radar  that  is  located 
throughout  the  United  States,  Alaska,  Hawaii,  and  Puerto  Rico.  The  WSR-88D  can  readily  detect 
birds  in  the  atmosphere  in  both  clear  air  and  precipitation  mode.  This  weather  radar  provides 
infonnation  on  movements  within  124  nautical  miles  (NM)  for  a  single  radar  station  as  well  as 
regional  and  national  scale  for  multiple  radar  sites  (Gauthreaux  and  Belser  1998).  Dokter  et  al 
showed  that  weather  radar  can  extract  near  real-time  bird  density  altitude  profiles  that  closely 
correspond  to  the  density  profiles  measured  by  dedicated  bird  radar  (Dokter,  et  al.  2011).  The 
WRS-88D  is  used  in  the  Avian  Hazard  Advisory  System  which  is  described  later  in  this  paper. 

Advisory  systems  provide  pilots  and  planners  a  method  for  estimating  risk  associated 
with  bird  strikes.  One  of  the  vanguards  of  these  advisory  systems  is  the  United  States  Bird 
Avoidance  Model  (USBAM).  The  USBAM  is  based  on  approximately  30  years  of  historic  bird 
observation  data  for  winter  and  summer  distributions.  These  point  data  were/are  transformed 
into  average  bird  mass  values  and  are  interpolated  spatially  in  a  GIS  environment  for  each  of  the 
bird  strike  relevant  species  with  a  resolution  of  1  km".  Between  the  winter  and  the  summer 
distribution  a  temporal  interpolation  is  conducted,  based  on  diurnal  and  annual  activity  pattern, 
breeding  success  and  mortality  rate.  The  overall  average  mass  of  birds  per  km  is  transformed 
into  bird  strike  risk  levels  of  low,  moderate  and  severe.  Model  output  is  displayed  in  an  internet 


map  application  that  combines  the  bird  strike  risk  level  infonnation  with  additional  important 
map  information  for  aviators. 

Modeling  in  this  case  stands  for  a  widely  automated  process  to  transform  the  historic  bird 
count  information  into  average  time  and  space  dependent  bird  strike  risk  levels.  Updates 
according  to  new  data  are  provided  approx,  every  2-5  years  (Ruhe  2005). 

The  risk  levels  that  USABAM  uses  describe  three  predicted  risk  classes  -  Low,  Moderate, 
and  Severe,  which  are  based  upon  the  bird  mass  in  ounces  per  square  kilometer.  The  "Moderate 
Zone"  indicates  a  risk  ratio  that  is  57-708  times  the  risk  of  the  "Low  Zone",  while  the  "Severe 
Zone"  indicates  a  risk  ratio  that  is  2,503-38,647  times  the  risk  of  the  "Low  Zone"  (United  States 
Air  Force  BASH  2012). 

The  species  data  was  acquired  from  several  key  datasets,  including  the  Audubon 
Societies'  Christmas  Bird  Count,  the  US  Biologic  Survey's  Breeding  Bird  Survey,  bird  refuge 
arrival  and  departure  data  for  the  contiguous  US,  and  many  additional  data  that  are  specific  to  a 
particular  bird  species. 

The  United  States  Avian  Hazard  Advisory  System  (AHAS)  is  an  online,  near  real-time, 
geographic  infonnation  system  (GIS)  used  for  bird  strike  risk  flight  planning  across  the 
continental  United  States.  Using  NEXRAD  (WSR-88D)  weather  radars  and  models  developed 
to  predict  bird  movement,  AHAS  monitors  bird  activity  and  forecasts  bird  strike  risk  as  well. 
AHAS  takes  current  weather  data  from  the  National  Weather  Service  and  calculates  the  risk 
large  bird  species  present,  based  upon  the  relationships  found  between  behavior,  weather  and 
strike  rate  with  each  species.  Standard  meteorological  calculations  are  used  to  detennine  thennal 
depth  and  strength.  Weather  data  is  also  used  to  detennine  when  birds  will  initiate  migration 
(Kelly,  et  al.  1999). 
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AHAS  consists  of  the  following  elements: 

•  A  forecast  of  bird  migratory  activity  for  the  next  twenty- four  hours. 

•  A  forecast  of  soaring  bird  activity  for  the  next  twenty-four  hours. 

•  Near  real  time  monitoring  of  bird  activity  with  NEXRAD  radar. 

•  Radar  Data  archiving  for  system  development. 

2.4.  Avian  Altitude  Density  Distributions 

For  bird  strikes  that  were  reported  at  <  500  feet,  passerines,  gulls  and  terns,  pigeons  and 
doves,  waterfowl,  and  birds  of  prey  were  the  species  groups  most  frequently  struck.  For  strikes 
above  500  feet,  waterfowl,  gulls  and  terns,  passerines,  birds  of  prey,  and  vultures  were  the 
species  groups  most  frequently  struck.  Waterfowl  comprised  53%  of  the  identified  birds  struck 
above  3,500  feet  (Dolbeer,  Height  Distribution  of  Birds  Recorded  by  Collisions  with  Civil 
Aircraft  2006).  Additionally,  74%  of  all  reported  strikes  were  <  500  feet  AGL. 

For  the  24%  of  strikes  above  500  feet,  and  using  height  interval  as  the  single 
independent  variable  and  number  of  reported  strikes  per  interval  as  the  dependent  variable, 
Dolbeer  determined  the  regression  equation  that  gave  the  best  correlation  for  reported  bird 
strikes  to  elevation  for  bird  strike  data  from  1990-2004.  His  analysis  resulted  in  a  negative 
exponential  model  with  height  as  the  independent  variable.  The  height  variable  explained  99% 
(R~  =  0.9891)  of  the  variation  in  number  of  bird  strikes  per  1,000-foot  interval  (starting  at  500 
feet). 

E[Strikes]  =  4469.2  exp(-0.3846  *  FlightAltitude)  ( 1) 

Plotting  this  function  against  altitude  gives  the  following: 
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Figure  1  Bird  Strike  Distribution 

Thus,  the  majority  of  strikes  happen  below  3500  feet  AGL.  However,  this  data  does  not 
explicitly  give  the  distribution  of  avian  altitudes. 

To  better  understand  avian  altitude  stratification,  specific  research  for  an  area  is  required. 
Bird  populations  are  constantly  changing  in  response  to  various  anthropogenic  as  well  as  natural 
factors.  It  is  estimated  that  the  longevity  of  the  predictions  by  bird  avoidance  models  (BAMs)  to 
be  in  the  order  of  5-10  yr.  Therefore,  a  model  update  is  recommended  approximately  once  every 
5  years  (Shamoun-Baranes,  et  al.  2008). 

In  the  Netherlands,  a  multidisciplinary  team  developing  a  Bird  Avoidance  Model  for 
Northwest  Europe  captured  data  on  60  different  species  of  birds  in  order  to  model  bird 
distribution  and  flight  altitude  predictors  and  distributions  (Shamoun-Baranes,  Sierdsema,  et  al. 
2005).  Data  was  collected  via  Hollandse  Signaal  Apparaten  (HSA;  Hengelo,  the  Netherlands) 
midlife  update  (MLU)-Flycatcher  tracking  radar.  Identification  and  further  refinement  was 
collected  using  two  methods:  1)  visual  identification  using  a  video  camera  with  a  300-mm  lens 
mounted  parallel  to  the  tracking  radar,  and  2)  the  classification  of  wing  beat  frequencies 
(Shamoun-Baranes,  Van  Loon,  et  al.  2006). 
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Multiple  linear  regression  models  were  built  by  fitting  explanatory  variables  to  maximum  hourly 
flight  altitude  as  the  response  variable.  Explanatory  variables  included  several  meteorological 
variables  such  as  temperature,  wind  speed  and  direction,  boundary  layer  height,  relative 
humidity,  lifted  index  (a  measure  of  atmospheric  instability),  time  of  day  and  time  of  year.  Each 
species  was  analyzed  separately.  In  order  to  keep  models  relatively  simple,  their  models  were 
limited  to  no  more  than  four  significant  (p  <  0.05)  explanatory  variables. 

The  following  represents  the  estimated  maximum  altitude  of  a  buzzard  based  upon  the 
explanatory  variables  of  relative  humidity  (RH),  boundary  layer  height  (BLH),  maximum  daily 
temperature  (Tmax),  and  lift  index  (LFTX). 

In  (height)  =  {-0.02[RH](%)  +  0.0001[BLH]  +  0.0n[Tmax]-0.04\{LFTX]  +  5.33}  (2) 

Additionally,  the  study  captured  the  altitude  distribution  of  the  species  studied.  Figure  2 
shows  the  distribution  for  the  buzzard  with  an  associated  gamma  distribution. 


Figure  2  Buzzard  Altitude  Distribution 

Multiple  studies  have  estimated  the  best-fit  parameters  for  the  observed  flight  height 
distribution.  The  empirical  data  was  fitted  using  maximum-likelihood  parameter  estimations. 
Various  literatures  (Intachat  and  Holloway  2000),  (Stumpf,  et  al.  2011)  and  (Shamoun-Baranes, 
Van  Loon,  et  al.  2006)  reported  that  the  gamma  distribution  (among  others)  was  an  appropriate 
distribution.  In  probability  theory  and  statistics,  the  gamma  distribution  is  a  two-parameter 
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family  of  continuous  probability  distributions,  with  a  shape  parameter  “a”  and  a  scale  parameter 
“(3”  called  the  rate  parameter.’ 

The  model  in  this  paper  uses  an  independent  gamma  distribution  (if  known)  for  each 
species  (if  known).  However,  if  the  species  of  bird  (but  not  the  percentages  of  each  species)  in 
an  AO  is  known  and  a  gamma  distribution  of  the  species  is  known,  then  one  must  explore  the 
ability  to  add  gamma  distribution  in  order  to  have  a  gamma  distribution  across  the  entire  AO.  If 
the  beta  parameters  of  the  gamma  distributions  are  equal,  then  the  gamma  distributions  can  be 
added. 

If  Xi,  X2,...,  Xn  are  independent  random  variates  with  a  gamma  distribution  having 
parameters  (ai,P),  (a2,P),...,(an,P)  then  £f=  1Xi  is  distributed  as  gamma  with  parameters: 

a  =  Yjai,P  =  p  (3) 

1=1 

If  however,  the  beta  parameters  are  different,  then  the  process  is  not  as  simple  as  outlined 
above.  In  this  situation,  it  is  feasible  to  1)  use  an  equal  percentage  for  the  known  birds  in  the 
AO  and  then  apply  the  specific  gamma  distribution  for  that  species  of  bird  or  2)  use  an 
approximation  of  the  parameters.  The  approximation  is  outlined  in  Thom’s  work  fitting  the 
gamma  distribution  for  precipitation  data.  Since  the  shape  (a)  parameters  add  in  convolution  of 
gamma  distributions,  and  the  mean  x  is  also  a  sufficient  estimator  of  p  (the  1st  moment)  it  is 
straightforward  to  determine  the  scale  ((1)  from  the  summed  means  (Thom  1968).  This  gives  the 
following  estimates  “g”  for  shape  and  “b”  for  scale: 

III 

g  =  YjSn  (4) 

1=1 

m  '  — 

b  -  V  — -  where  m  -  #  of  gamma  distributions 

/=1  & 
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2.5.  Independence  -  Bird  Flocks  /  Small  Bird  Estimates 


Because  radar  is  primarily  utilized  to  detect  and  identify  bird  species  in  an  area  of 
operations  common  to  low  level  routes  or  bombing  ranges,  the  flocking  behavior  of  birds  can 
make  it  difficult  to  determine  the  number  of  birds  in  the  flock  and  thus  the  type  of  bird  present. 
One  technique  that  many  researchers  utilize  is  the  “small  bird  estimate”  that  was  first  utilized  by 
Kelly  (1995).  Birds  of  each  size  are  scaled  to  standardize  the  birds  by  mass  into  categories.  The 
use  of  small  bird  estimates  (SBE)  helps  to  counter  the  problem  of  unknown  bird  numbers  per 
radar  target.  A  medium  target  on  the  radar  screen  may  be  a  single  intermediate-sized  bird  or  a 
small  flock  of  small  birds.  Either  way  it  is  represented  in  the  researcher’s  model  as  the  same 
number  of  SBE’s.  Their  assumption  is  that  it  is  equally  hazardous  to  strike  one  intermediate¬ 
sized  bird  or  a  small  flock  of  small-sized  birds.  A  larger  flock  of  small  birds,  a  small  flock  of 
intermediate-sized  birds,  and  an  individual  large  bird  would  all  be  categorized  as  large  bird 
targets  and  would  be  recorded  as  the  same  number  of  SBEs.  Thus,  the  numbers  birds  per  radar 
target,  hence  risk,  though  not  completely  quantifiable,  is  incorporated  in  the  algorithm 
(Zakrajsek  and  Bissonette  2002). 

Budgey’s  research  indicated  that  there  is  an  apparent  relationship  between  the  wingspan  of  a 
flocking  bird,  its  nearest  neighbor  distance  and  the  number  of  birds  likely  to  be  struck  by  an 
aircraft  encountering  a  flock  of  that  species  (Budgey  1998).  His  research  was  limited  in  that  not 
all  species  were  studied  and  the  number  of  flocks  per  species  studied  was  small.  If  his  research 
holds  for  all  species,  then  this  method  of  determining  the  number  of  birds  in  a  flock  and  the 
probability  of  the  number  of  bird  strikes  could  be  used  to  enhance  current  collision  models.  One 
would  need  to  detennine  the  probability  of  encountering  the  flock,  and  then  determine  the 
probability  of  a  strike  from  within  the  flock. 
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For  the  models  used  in  this  paper,  it  will  be  assumed  that  the  bird  species  identification  is  known 
either  through  visual,  historical  or  SBE  approximation.  Flocking  behavior  and  numbers  per 
flock  will  not  be  utilized.  If  there  is  enough  fidelity  in  detennining  birds  in  an  AO,  then  each 
bird  will  be  treated  independently  and  assumed  to  follow  a  Poisson  Spatial  distribution. 

Risk  Level  Ranking 

All  wildlife  species  are  not  equally  hazardous  to  aviation.  The  hazard  of  a  species  to  an 
aircraft  is  strongly  related  to  the  species  mean  body  mass.  However,  it  is  not  the  only  factor 
when  modeling  the  hazard  associated  with  a  specific  bird. 

The  FAA  historical  model  looked  at  the  hazard  level  for  108  bird  species  with  25  or  more 
strikes,  based  on  the  percentage  of  strikes  that  have  caused  an  adverse  effect  (damage  and/or 
negative  effect  on  flight),  ranged  from  0  percent  for  8  species  to  80  percent  for  snow  geese  (FAA 
2012). 

The  Royal  Netherlands  Air  Force  models  their  bird  hazard  index  to  be  equal  to  the 
relative  strike  sensitivity  times  the  damage  probability. 

(Morgenroth  2003)  modeled  bird  strike  hazards  based  on  behavioral  specific  aspects  of 
birds.  Zakrajsek  and  Bissonette’s  (2002)  model  was  based  on  military  bird  strike  statistics  and 
specifically  focused  on  the  number  of  damaging  strikes  and  the  cost  associated  with  that  damage 
as  a  criteria. 

DeVault  used  percentage  of  strikes  with  damage,  percentage  of  strikes  with  substantial 
damage  and  percentage  of  strikes  with  effect  on  flight  (EOF)  as  factors  in  which  species  were 
ranked  and  a  relative  hazard  score  was  calculated.  For  birds,  they  assessed  effects  of  body  mass, 
body  density,  and  group  size  on  relative  hazard  scores  for  77  species.  Additionally,  they  only 
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used  data  from  strikes  that  occurred  at  less  than  500  feet  AGL  (DeVault,  et  al.  2011).  Table  2 


highlights  the  top  15  hazardous  birds  based  upon  their  research  and  model. 


Table  2  DeVault  Top  15  Relative  Hazard  Score 


Species 

Tota  1 

strikes 

reported 

Composite 

ra  nk 

Relative 

hazard 

score 

Body  mass 

(g) 

%  of  strikes  with 

mult,  bi rds 

Other  geese* 

20 

1 

100 

2290 

60 

Canada  goose 

776 

2 

76 

3564 

47.9 

Other  ducks* 

77 

2 

78 

916 

46.8 

Turkey  vulture 

159 

2 

73 

1467 

9 

Double-crested 

24 

5 

71 

1674 

16.7 

Great  horned  owl 

29 

5 

72 

1309 

3.4 

Brown  pelican 

31 

7 

66 

3348 

9.7 

Sandhill  crane 

66 

8 

61 

5571 

44.6 

Glaucous-winged  gull 

27 

9 

64 

1010 

25.9 

Wild  turkey 

38 

9 

65 

5811 

23.7 

Bald  eagle 

74 

11 

59 

4740 

12.2 

Great  black-backed  gull 

20 

12 

53 

1659 

15 

Osprey 

77 

13 

53 

1485 

2.6 

Great  blue  heron 

132 

14 

51 

2390 

2.3 

Ring-necked  pheasant 

45 

15 

47 

1135 

8.9 

Dolbeer,  Wright  and  Cleary  created  a  model  to  rank  19  different  bird  species  based  upon 
1)  Damage,  2)  Major  Damage,  and  3)  Effect  on  Flight.  Using  rankings  within  each  category, 
they  developed  a  composite  ranking  and  a  relative  hazard  score.  The  relative  rank  of  each 
species  group  was  compared  with  every  other  group  for  the  3  categories,  placing  the  species 
group  with  the  greatest  hazard  rank  for  >  2  of  the  3  categories  above  the  next  greatest-ranked 
group.  Their  relative  hazard  score  was  related  strongly  to  mean  body  mass  for  the  19  bird 
species  groups  (R~  =  0.71,  p<.01).  Table  3  contains  their  final  results  (Dolbeer,  Wright  and 
Cleary,  Ranking  the  Hazard  Level  of  Wilflife  Species  to  Aviation  2000). 
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Table  3  Risk  Level  Ranking 


Ranking  by  Criteria 

Species  Group 

Damage 

Major 

Da  mage 

Effect  on 
flight 

Composite 

Ranking 

Relative 

Haza rd  Score 

Deer 

1 

1 

1 

1 

100 

Vultures 

2 

2 

2 

2 

63 

Geese 

3 

3 

4 

3 

52 

Cra  nes 

4 

4 

7 

4 

48 

Osprey 

6 

5 

3 

5 

50 

Pelica  ns 

5 

7 

5 

6 

44 

Ducks 

7 

6 

8 

7 

37 

Hawks 

8 

15 

9 

9 

25 

Eagles 

8 

15 

9 

9 

31 

Rock  Dove 

11 

8 

11 

10 

24 

Gulls 

10 

11 

13 

11 

22 

Herons 

12 

14 

12 

12 

22 

Mourning  Dove 

14 

9 

17 

13 

17 

Owls 

13 

12 

19 

14 

16 

Coyote 

15 

17 

6 

15 

20 

American  Kestrel 

16 

10 

16 

16 

14 

Shorebi  rds 

17 

19 

14 

17 

12 

Crows  -  Ravens 

18 

16 

15 

18 

12 

Blackbirds  -  Starling 

19 

18 

18 

19 

9 

Sparrows 

20 

21 

20 

20 

4 

Swallows 

21 

20 

21 

21 

2 

2.6.  Spatial  Poisson  Process 

Spatial  point  pattern  data  occurs  frequently  in  a  wide  variety  of  scientific  disciplines, 
including  seismology,  ecology,  forestry,  geography,  spatial  epidemiology,  and  material  science 
(Moller  and  Waagepeterson  2006).  One  may  think  of  a  spatial  point  process  as  a  random 
countable  subset  of  a  space  X.  Assume  that  X  c  Typically,  X  will  be  a  d-dimensional  box 
or  all  ofRd.  However,  it  could  also  be  Xd_1,  the  (d-1)  -  dimensional  unit  sphere  (Johnson  2010). 
In  many  instances,  one  may  only  observe  points  in  a  bounded  subset  (window)  A  Q  X.  For 
example,  there  may  be  interest  in  the  spatial  distribution  of  a  certain  species  of  tree  in  a 
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wilderness  area.  Without  having  to  count  or  record  every  tree,  efforts  should  be  concentrated  on 
several  square  windows  where  A;  =  [a,o,  an]  ,  an  >  aio,^4j  fl  Aj  =  0,  j  A  i  and  A  =  At  c  X. 
One  of  the  simplest  and  fundamental  spatial  point  processes  is  the  completely  random,  or  spatial 
Poisson  point  process  (Guttorp,  Brillinger  and  Schoenberg  2002).  Spatial  Poisson  processes  play 
a  fundamental  role  in  the  theory  of  point  processes.  They  possess  the  property  of  “no 
interaction"  between  points  or  “complete  spatial  randomness".  The  points  are  stochastically 
independent  and  the  probability  of  the  number  of  points  N(A)  in  a  region  A,  is  given  by  the 
Poisson  distribution  as  outlined  below: 

Let  QA  —  {^-(x,.)^  <z  X,  N  e  Aujoo}}  denote  the  space  of  configurations  on 

x  =  ad,d>  l. 


Figure  3  Spatial  Poisson  Process  in  2-Space 
Let  ®(d)  =  #{xe©:xed}  =  H.JP')- 

The  Poisson  Probability  measure  Pa  with  intensity  p(x)dx  on  X  satisfies 

v  ...  (o|A|)n 

Hj  (go  G  nx  :  to(A)  =  n)  =  e  - - — ,  n  6  hi 

with  a\A\  =  jAp(x)dA  =  fRdlA(x)p(x)dx. 

Usually  the  intensity  function  p(x)  will  be  constant  (homogeneous),  i.e.  p(x)  =  1  >  0,  x  £  X, 
where  X  >  0  is  called  the  intensity  parameter: 
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A  practical  application  of  the  spatial  Poisson  process  (SPP)  distribution  is  observed  in  a 
minefield  crossing  (Kim  2002).  In  this  example,  the  premise  is  that  mines  are  buried  in  a  field  of 
known  dimension.  The  location  of  the  mines  is  not  known,  but  the  number  of  mines  in  this 
bounded  area  is  known.  The  author  asserts  that  the  placement  of  the  mines  follows  a  SPP  and 
they  will  detonate  if  you  are  with  a  certain  distance  “r  =1.5  feet”  from  the  mine.  Suppose  the 
field  is  100  ft  x  200  ft  and  you  are  crossing  it  at  its  shortest  distance  (100  ft).  The  intensity  of  the 
mines  is  X  =  15/20000  =  .00075  per  square  foot. 

not  to  scale 

100  ft 

200  ft  ' — r — J 

2*r  =  3  feet 

Figure  4  Minefield  Example 

The  region  |A|  in  this  case  is  3  ft  x  100  ft  =  300  square  ft.  The  probability  of  safely  transiting  the 
minefield  (probability  of  no  strikes)  is: 

p0  =  P(o>  g  :  oj(A)  =  0)  =  (0-0007^*300)  g-o. 00075*300  =  0J985 16  (6) 
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Thus,  in  this  situation,  there  is  a  -80%  chance  of  traversing  the  minefield  safely. 

This  concept  can  easily  be  extended  to  3-space.  In  the  following  example,  the  probability  of  an 
aircraft  /  bird  encounter  (defined  as  a  bird  being  within  a  certain  proximity  to  the  aircraft)  is 
desired.  Table  4  contains  bird  sightings  from  Dover  AFB  (KDOV)  during  the  morning  of  2 
October,  2011. 


Table  4  Dover  AFB  Bird  Sightings 


Date 

Ti  me 

Zone 

Species  ID 

Number 

Observed 

Habitat 

Behavior 

Observation 

Method 

10/2/2011 

08:29:31 

Ai  rfi  e  1  d 

Canada 

Geese 

16 

Sky 

Overflight 

Visua  1 

10/2/2011 

08:27:49 

Ai  rfi  e  1  d 

Gulls 

80 

Runway 

Loafi  ng 

Visua  1 

10/2/2011 

08:38:49 

Ai  rfi  e  1  d 

Gulls 

62 

Runway 

Loafi  ng 

Visua  1 

10/2/2011 

08:48:50 

Ai  rfi  e  1  d 

Gulls 

4 

Sky 

Overflight 

Visua  1 

In  this  scenario,  a  C-5  Galaxy  departing  Dover  AFB  is  modeled  with  the  assumption  that  the 
birds  are  distributed  in  the  airfield  environment  according  to  a  SPP  and  the  speed  of  the  birds  is 
negligible  compared  to  the  speed  of  the  aircraft  (thus  they  are  considered  stationary).  The 
encounter  space  |A|  is  the  volume  of  airspace  around  the  aircraft  equal  to  the  radius  of  the 
wingspan  swept  through  the  distance  the  aircraft  flies  in  the  AO.  This  results  in  a  disk  of  radius 
=  1 1 1  ft  being  swept  through  the  air  for  the  distance  of  the  takeoff  roll  and  for  the  distance 
required  to  climb  out  of  the  airport  environment  (in  this  case  500  ft  of  altitude).  Table  5  contains 
the  data  for  this  scenario. 
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Airfield  Environment  Volume 

"AO" 

Circle  with  diameterof  12903  ft  (length  of  runway) 

Altitude  =  500  ft  (approx  altitude  when  passing  departure  end  of 

runway) 

( p i ) ( 6451.5) ( 645 1.5) ( 500) =65,346,308,032 .5  ft  cubed 

C-5  Wingspan 

222  ft 

Takeoff  Roll 

8400  ft 

Takeoff  Roll  "Volume" 

(pi )(111)(111)(8400)=324,948,696  ft  cubed 

Climb  Rate  /  Climb  Distance 

1800  ft/sec/4527  ft 

Climb  "Volume" 

(pi )(111)(111)(4527)=175, 140,304  ft  cubed 

X 

162/65346308033  =  2.4791E-09 

1 A  | 

324,978,696  +  175,140,304.4  =  500119000  ft  cubed 

Table  5  SPP  in  3  Space  Example  Data 


Figure  5  is  a  graphical  depiction  of  the  scenario: 

6451.5  ft 


Figure  5  SPP  in  3  Space  Depiction 


Using  Equation  5  with  the  data  from  Table  5,  the  resulting  probability  of  no  (zero)  aircraft  /  bird 
encounters  is  -28.94%.  The  probability  of  one  (1)  encounter  is  -35.88%.  The  probability  of  one 
or  more  strikes  is  1  -  0.2894  =  0.7106  or  -  71%.  By  using  the  disc  of  radius  1 1 1  ft  (the 
wingspan),  this  model  overestimated  the  area  being  swept  through  the  AO. 

This  example  assumed  that  the  relative  velocity  between  bird  and  aircraft  as  negligible 
and  it  assumed  that  the  birds  were  stationary  during  the  encounter.  A  commonly  held  belief 
(aircrew  legend)  is  that  birds  will  dive  if  they  encounter  an  aircraft.  However,  research  has 
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shown  that  the  most  important  feature  of  the  execution  of  last-minute  avoidance  maneuvers  by 
birds  has  been  shown  to  be  their  unpredictability.  Avoidance  of  aircraft  in  slow,  straight  and 
level  flight  seems  to  present  little  difficulty  but  at  the  speed  and  low  altitudes  at  which  birds  are 
most  likely  to  be  encountered,  aircraft  are  usually  climbing  or  descending.  This  is  where  the 
unpredictability  comes  in.  Observational  evidence  suggests  that  the  tendency  to  attempt  to  dive 
or  free  fall  beneath  an  aircraft  rather  than  climb  above  it,  which  is  more  marked  for  some  species 
than  others,  is  not  reliable  enough  to  support  a  corresponding  avoidance-climb  response.  Such  a 
response  might  still  lead  to  a  strike,  while  the  aircraft  would  be  either  at  an  increased  rate  of 
climb  or  in  the  transition  from  descent  to  climb,  and  aircraft  performance  margins  would  be 
reduced  with  a  greater  risk  of  significant  consequences  arising  from  strikes.  Part  of  the 
underlying  explanation  for  unpredictable  aircraft  avoidance  is  that  birds  do  not  always  seem  able 
to  perceive  aircraft  as  being  in  motion.  Exceptionally,  some  hawks  and  eagles  have  been  seen  to 
‘attack’  an  aircraft  which  they  encounter  (SKYbrary  2012).  Therefore,  continue  to  assume  the 
stationary  location  of  birds  with  respect  to  aircraft  for  modeling. 

2.7.  Related  Random  Collision  Models  for  Avians 

Several  researchers  have  studied  the  expected  numbers  of  individual  prey  captured  (#  of 
birds  ingested  into  an  aircraft  engine)  from  flocks  of  birds  (Major,  Dill  and  Eaves  1986).  Their 
study  determined  the  expected  number  of  collisions  for  various  predator  speeds  and  trajectories, 
flock-predator  initial  distances  and  angles,  and  flock  sizes,  shapes,  densities,  trajectories,  and 
speeds.  Generally,  larger  predators  and  clustered  predators  (bigger  engines  in  a  triangular  vice 
straight  line  configuration)  caught  more  prey.  A  moving  three-dimensional  Poisson  model  was 
chosen  to  describe  the  interactions  between  predators  (sets  of  engine  frontal  areas)  and  prey 
(dunlin  flocks).  The  model,  incorporated  into  an  interactive  computer  program,  calculated  the 
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expected  number  of  captures  of  birds  by  each  predator  for  each  operator-supplied  set  of  predator- 
prey  trajectories  and  speeds.  The  model  was  also  used  to  calculate,  for  individual  predators,  the 
exponential  distribution  of  the  inter-bird  capture  interval,  its  expected  value,  and  its  variance 
(Major,  Dill  and  Eaves  1986). 

Another  area  of  study  with  respect  to  random  collisions  of  birds  and  objects  in  motion  is 
wind  turbine  rotors.  When  a  bird  flies  through  the  disc  swept  out  by  blades  of  a  wind  turbine 
rotor,  the  probability  of  collision  depends  on  the  motions  and  dimensions  of  the  bird  and  the 
blades.  V.A.  Tucker  created  a  mathematical  model  that  predicts  the  probability  of  a  collision 
when  a  bird  flies  through  the  disk  swept  out  by  the  rotor  blades  (Tucker  1996). 

2.8.  Extended  Spatial  Poisson  Process 

Extended  Poisson  process  models  derive  their  name  from  the  fact  that  they  are  based 
upon  generalizing  the  simplest  Markov  birth  process,  the  Poisson  process  (T.  Holzmann  2009). 
They  involve  representing  a  discrete  distribution  as  the  distribution  of  the  number  of  events 
occurring  in  a  finite  time  interval  of  a  state-dependent  Markov  birth-death  process  (Podlich, 
Faddy  and  Smyth  1999)  &  (Faddy  1997). 

Holzmann  and  Cochran  furthered  Faddy’s  (1990  &  2008)  study  of  spatial  data  with 
respect  to  Extended  SPP.  They  utilized  a  linearly  decreasing  arrival  rate  for  the  birth-death 
process  of  the  Markov  transition  in  order  to  “remove”  encounter-able  elements  from  the  space 

(Holzmann  and  Cochran  2012).  This  resulted  in  A'  —  —(N  —  n)ln  (1  —  p)  where  p  is  the 

1 

density  function  (calculated  by  —  ,  where  |V|  is  the  “size”  of  the  AO)  and  N  is  the  number  of 
encounter-able  entities  in  the  AO.  This  leads  to  the  Markov  transition  matrix  Q  given  by: 


23 


-N; l  NX  0  •••  0 

0  -(A-l)/l  (iV-l)2  0 

Q=  0  0  -(N-2)X  0  (7) 

0  0  0  0  0 

where  A  =  In  (1  — 

This  then  leads  to  the  probability  distribution  for  the  number  of  encounters  in  the  AO  for  the 
specific  window  |A|  as: 

V  =  PoexP  (<2  *  \A\) 


Where:  p0  =  {1,  0, 0,  ,  0]  is  the  initial  condition  of  no  (0)  encounters  and  the  probability  pt  of 

encounters  at  a  count  size  of  i  (i  =  0,  1,. . .,N)  is  the  (/  +  l)th  element  of p. 


Example: 

Given  a  space  W  of  size  |V|  with  5  entities  distributed  throughout  it,  an  object  A  passes 
through  space  W  sweeping  a  path  of  size  |A|. 


|V|  =  10000  units  cubed 
|A|  =  1500  units  cubed 
N  =  5 

The  following  is  the  transition  matrix: 

0.095  0.284  0.346  0.205  0.063  0.007 

0  0.152  0.366  0.330  0.132  0.020 

0  0  0.243  0.439  0.264  0.053 

0  0  0  0.390  0.469  0.141 

0  0  0  0  0.624  0.376 

0  0  0  0  0  1.000 
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And 


p  =  [0.095  0.284  0.346  0.205  0.063  0.007] 

Therefore,  the  probability  of  zero  (0)  encounters  of  object  A  with  an  entity  is  -9.5%,  and  the 


probability  of  one  (1)  encounter  is  -28.4% 

However,  if  one  were  able  to  detect  and  “remove”  an  entity  from  the  encounter-able  list, 
then  the  initial  condition  could  be  set  at  p±  =  {0, 1,  0,  0, 0, 0]  and  row  2  of  Equation  8  would  be 
used.  The  Extended  Spatial  Poisson  Process  has  application  in  the  avian  encounter  models  when 
calculating  a  risk  score  because  it  can  “remove”  an  entity  that  has  already  been  encountered  and 
therefore  can  eliminate  double  counting  and  have  a  more  realistic  model  by  assuming  that  a  hit 
bird  is  removed. 


2.9.  Matrix  Exponentiation 

The  calculation  of  the  matrix  exponential,  or  an  approximation  of  p,  is  crucial  to  the 
Extended  Spatial  Poisson  Process.  However,  calculating  the  matrix  exponential  is  not 
straightforward.  The  most  commonly  used  algorithm  to  calculate  the  exponential  is  that  of 
Golub  and  Van  Loan  (Faddy  and  Smith,  Extended  Poisson  process  modeling  of  dilution  series 
data  2008).  Current  applications,  such  as  ®MATLAB  and  ®SAS  quickly  solve  matrix 
exponentiation  of  reasonable  size.  However,  these  programs  are  not  readily  available.  Appendix 
II  contains  VBA  code  that  allows  the  user  to  obtain  p  by  an  iterative  method  of  Taylor 
expansions: 


(9) 
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Additionally,  it  also  worth  noting  that  Q  *  |^|  (referred  to  a  Q  from  here  on  out)  has  no 
repeated  eigenvalues  (the  diagonal  elements  are  unique)  and  Q  is  upper  triangular.  Given  this 
information,  Q  can  be  definitely  be  diagonalized. 

If  the  eivenvalues  of  Q  are  Xi,  X2,...,  Xn,  then  each  eigenvector  (xi,  X2,...,  xn)  of  Q  is 
linearly  independent.  If  these  eigenvalues  are  the  columns  of  matrix  S,  then  S'1  Q  S  is  the 
diagonal  matrix  A  and  the  eigenvalues  of  Q  are  the  diagonal  elements  of  A.  This  result  in  the 
following  diagonal  matrix:  (Strang  2006). 


Q  =  SAS  1  ^SlQS  =  A 


K 


A. 


Additionally,  it  follows  that: 

eQ=S-eQ-S~l=S 


(10) 


(11) 


Therefore,  the  eigenvectors  associated  with  the  eigenvalues  of  Q  can  be  determined  and  Q  can  be 
exponentiated  as  well. 
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3.  Methodology 


3.1.  Overview 

This  chapter  describes  the  formulation  and  use  of  the  model.  The  Spatial  Poisson 
Process,  using  specified  parameters,  will  be  used  to  determine  the  probability  of  an  aircraft  /  bird 
encounter.  This  probability  will  then  be  aggregated  within  the  specific  bird  species  and  then 
multiplied  by  a  scaled  risk  score  and  then  all  species  will  be  summed  in  order  to  give  the 
decision  maker  an  additional  analysis  tool  when  planning  missions. 

3.2.  Assumptions 

The  following  assumptions  are  made  with  respect  to  model  formulation  and  use: 

•  The  baseline  for  an  Operations  Area  (AO)  is  1km  x  1km  x  3500  feet 

•  Estimated  number  and  type  of  birds  in  local  area  are  known  to  an  acceptable  level 

•  The  velocity  of  birds  in  the  AO  is  negligible  compared  to  the  speed  of  the  aircraft  and  can 
be  considered  stationary 

•  The  number  of  birds  in  the  AO  is  constant 

•  An  aircraft/bird  encounter  within  the  ellipse  formed  by  the  height  and  width  of  an  aircraft 
will  be  considered  a  strike 

•  Altitude  distributions  of  birds  (if  known)  will  following  a  gamma  distribution  (birds 
within  an  altitude  band  will  be  distributed  via  a  Poisson  process) 

•  SBE’s,  as  outlined  in  the  independence  section,  are  utilized  as  needed  to  detennine 
species  and  numbers  of  birds  in  the  AO. 

3.3.  Methods 

The  Poisson  spatial  srocess  in  3-space  will  be  utilized  to  determine  the  probability  of  a 
bird  strike.  Recall  Equation  5: 
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d/  r\x  t  /<\  \  (^M|) 

pn=P(coeQ.  :o)(A)-n)- — — — e  1 

n\ 

|A|  is  the  volume  of  space  defined  by  sweeping  the  elliptical  disc  around  an  aircraft 
through  the  AO  for  a  defined  distance  d.  The  ellipse  around  an  aircraft  is  a  greater  area  than  the 
frontal  area  of  an  aircraft  and  can  be  substituted  for  actual  frontal  area  in  order  to  provide  a  more 
precise  calculation.  The  following  is  an  illustration  of  the  ellipse. 


k  is  the  intensity  function  =  #  birds  in  the  AO  /  Total  volume  of  AO 
Lambda  will  be  scaled  by  the  following  parameters: 

•  gamma  height  distribution  of  the  specific  bird  species 

•  the  percentage  of  a  specific  bird  species  to  the  total  number  in  the  AO. 

•  “n”  is  the  number  of  birds  in  the  AO 

The  Dolbeer  risk  score  ranking  (Table  2)  will  be  used  for  the  risk  each  species  presents. 
This  data  is  scaled  using  Equation  12  (removing  deer  and  coyote)  and  the  result  is  the  multiplier 
for  which  the  “probability  of  bird”  strike  will  be  multiplied  against  in  order  to  determine  a  risk 
score. 
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Species 

Relative 

Haza rd  Score 

Vultures 

63 

Geese 

52 

Cra  nes 

48 

Osprey 

50 

Pelicans 

44 

Ducks 

37 

Hawks 

25 

Eagles 

31 

Rock  Dove 

24 

Gulls 

22 

Herons 

22 

Mourni ng  Dove 

17 

Owls 

16 

American  Kestrel 

14 

Shorebi  rds 

12 

Crows  -  Ravens 

12 

Blackbirds  -  Starling 

9 

Sparrows 

4 

Swa  1  lows 

2 

Figure  6  Risk  Score  Multiplier 

The  methods  used  and  the  parameters  needed  result  in  three  cases  for  the  planners.  The 
baseline  inputs  for  all  three  cases  consist  of  complete  AO  data  and  aircraft  data.  The  following 
are  the  specific  case  inputs  and  outputs. 

Case  #1  inputs  are  total  birds  in  the  AO.  Case  #1  outputs  are  the  probabilities  of  an 
encounter. 

Case  #2  inputs  are  total  birds  in  the  AO  and  the  specific  species  percentages  in  the  AO. 
This  results  in  outputs  of  the  probabilities  of  an  encounter  and  a  risk  score  /  risk  matrix. 

Case  #3  inputs  are  total  birds  in  the  AO  and  the  specific  species  percentages  in  the  AO, 
gamma  height  distribution  and  the  aircraft  altitude  band  to  be  flown.  The  outputs  for  case  #3  are 
the  probabilities  of  an  encounter  and  a  risk  score  /  risk  matrix,  with  a  higher  fidelity  due  to  the 
addition  of  the  gamma  distribution  and  aircraft  altitude  band. 
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3.4.  Parameters  and  Calculation 


The  following  are  user  inputs: 


Identifier 

Description 

Used  in 

L 

AO  Length  in  meters 

V  A  AO  Volume 

W 

AO  Width  in  meters 

V  A  AO  Volume 

H 

AO  Height  in  meters 

V  A  AO  Volume 

A_Mx 

Aircraft  Altitude  Max 

v  A  Modified  AO  Volume 

A_Mn 

Aircraft  Altitude  Min 

v  A  Modified  AO  Volume 

A_W 

Aircraft  Wingspan 

A  A  Aircraft  Volume 

AH 

Aircraft  Height 

A  A  Aircraft  Volume 

Surface  Area 

Aircraft  Frontal  SA 

A  A  Aircraft  Volume 

AD 

Aircraft  Distance 

A  A  Aircraft  Volume 

T 

Total  #  birds  in  AO 

X  A  Intensity  Function 

®  species^1 

Gamma  Alpha  value 

L’spedesAlt  A  Intensity  Function  modified  volume 

-^speciesP 

Gamma  Beta  value 

L’spedesAlt  A  Intensity  Function  modified  volume 

R  T 

0  species 1 

%  of  T 

A’spcccs^  Intensity  Function  bird  %  known 

R  RT 

^“species 

Species  Risk  Score 

Rspecies  —  total  risk  by  species 

The  following  calculations  are  needed  for  model  implementation: 


|V|  =  L  *  W  *  H 
v  =  L  *  W  *  (A_Mx  -  A_Mn) 

|A|  =  1/2(A_W)  *  14(A_H)*  n  *(A_D) 


{length  *  width  *  height} 

{length  *  width  *  altitude  band  height} 
{area  of  ellipse  *  distance} 


or 


|A|’  =  Surface_Area*(A_D) 
A,  =  T  /  |V| 

=  (B  species  T  *  T  )  /  |V| 


Vs 

V 


{Aircraft  Frontal  Surface  Area  *  distance} 


species  £ 


Alt  =  [(B specieS T  *  T  ) * (GammaCDF (A_Mx)  -  GammaCDF(AMn)]  /  |v| 
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3.5.  Flow  Diagram 


The  following  flow  diagram  illustrates  the  inputs,  decisions,  calculations  and  outputs  for 
the  model. 


Figure  7  Model  Flow  Chart 
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3.6.  Model 


This  section  describes  the  model  formulation.  There  are  three  different  models  that  can  be 
used  to  determine  a  bird  strike  and  an  associated  risk  score,  based  upon  the  number  and  type  of 
known  parameter  inputs.  These  models  use  the  spatial  Poisson  process  to  determine  a  strike 
probability.  The  strike  probability  is  used  to  detennine  the  expected  number  of  hits  for  a  flight 
through  an  AO  for  a  specific  species  of  bird.  The  expected  number  of  hits  (subtracting  double 
counting  /  intersections)  per  specific  species  is  then  multiplied  by  the  risk  multiplier  to  get  a  risk 
score.  These  individual  risk  scores  are  then  aggregated  over  all  the  species  to  get  a  total  risk 
score.  Figure  7  depicts  the  interaction  of  inputs,  calculations  and  outputs  for  the  models. 

The  change  in  the  input  parameters  affects  the  X  used  for  the  probability  function  and 
thus  the  strike  probability  and  the  risk  score.  The  following  table  outlines  the  inputs  and  outputs 
for  the  three  cases. 


Table  6  Model  input/output 


INPUT 

OUTPUT 

Case  # 

AO 

Volume 

Ai  rcra ft 

Distance 

Through 

AO 

Ai  rcra  ft 

Altitude 

Ba  nd 

Total  # 

Bi  rds  i  n 

AO 

Specific 
Species  % 

Avian 

Strike 

Probability 

Risk  Score 

1 

X 

X 

X 

X 

X 

2 

X 

X 

X 

X 

X 

X 

X 

3 

X 

X 

X 

X 

X 

X 

X 

X 

X 

The  following  equations  are  used  for  the  three  cases  outlined  in  Table  6  and  fonn  the  basis  for 
the  model.  These  equations  are  embedded  in  the  decision  support  system  in  order  to  generate 
strike  probability  and  risk  score. 

Case  #  1:  General  encounter  model: 

The  probability  of  an  encounter  for  all  the  birds  in  the  AO  is  calculated  using  the  spatial 
Poisson  process: 
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Pn  = 


g-A*|A| 


(A  *  |A|)n 
n! 


This  case  gives  the  probability  of  an  encounter  with  only  the  basic  parameters  of  AO  size, 
aircraft  size  (surface  area  or  wingspan/height),  aircraft  distance,  and  total  number  of  birds 
known.  There  is  no  risk  score  associated  with  this  case  because  specific  species  percentages  are 
not  known. 

Case  #2:  General  encounter  model  with  risk  score  added: 

The  probability  of  an  encounter  for  a  specific  species  is  calculated  using  the  Spatial 
Poisson  Process  with  a  modified  intensity  parameter  A'species  (generated  using  percentage  of 
specific  bird  in  the  AO). 

f  _y  .  *|A|  (^species  *  l^l) 

Pspecies(n)  =  e  |A|  * - - - 


The  risk  score  (by  species)  is  calculated  by  aggregating  over  the  specific  species  (the  number  of 
birds  x  probability  of  hit  (eliminating  double  counts))  x  species  specific  risk  level: 

N 


^species  —  ^  ’  tl  *  Pspecies(^)  *  RH 


species 


n=0 


Total  Risk  Score  is  calculated  by  summed  over  all  the  species. 

Rspecies 

All  species 


This  model  gives  the  probability  of  an  encounter  for  each  specific  species  and  a  risk  level 


score. 

Case  #  3:  General  encounter  model  with  risk  score  for  a  given  altitude  band: 

This  case  adds  another  element,  the  gamma  height  distribution,  to  the  problem  set.  By 
using  the  gamma  distribution  in  the  intensity  parameter  (A'species)  function,  the  intensity 
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parameter  now  becomes  itself  a  random  variable.  This  is  a  doubly-stochastic  process,  where  the 
intensity  parameter  (A'species)  is  now  also  stochastic  process.  This  is  also  called  a  Cox  process 
(Cox  1955). 

The  probability  of  encounter  within  the  altitude  band  is  detennined  using  a  doubly- 
stochastic  process  modified  intensity  parameter  A'speciesAlt: 

PspeciesAlt(n)  = 

The  risk  score  (by  species)  is  calculated  by  aggregating  over  the  specific  species  (the  number  of 
birds  x  probability  of  hit  (eliminating  double  counts))  x  species  specific  risk  level: 

N 

RspeciesAlt  —  I  tt  *  PspeciesAlt(n)  *  RHspecies 
n=0 


Total  Risk  Score  is  calculated  by  summing  over  all  the  species. 


R  = 


I 

All  species 


R 


species 


Alt 


This  model  calculates  the  probability  of  an  encounter  (eliminating  double  counts)  and  the 
risk  score  for  a  particular  band  of  altitude  within  the  overall  AO.  The  model  requires  a  gamma 
distribution  for  the  specific  species  that  are  in  the  AO  and  the  aircraft  altitude  band  to  be  flown  in 
order  to  calculate  these  values.  Figure  8  is  a  visual  representation  of  the  case  #3. 
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Figure  8  Gamma  Altitude  Distribution  and  Altitude  Band 
In  order  to  give  a  more  precise  estimate  and  to  eliminate  entities  that  have  already  been 

encountered,  one  can  apply  the  extended  spatial  Poisson  process  to  the  encounter  probability 
when  determining  the  risk  score.  Therefore,  the  model  detennines: 

P(x+1  |  x  ,  where  x  =  1,  2,  ...,  N-l  &  N  =  total  #  birds  in  AO). 

Figure  9  and  Figure  10  are  the  outputs  from  the  spatial  Poisson  process  and  the  extended 
spatial  Poisson  process  respectively  for  the  same  scenario.  In  this  scenario,  the  risk  score  =  1  for 
brevity. 


n 

P(n) 

Risk 

0 

0.675232 

0 

1 

0.265163 

0.26516 

2 

0.052065 

0.07652 

3 

0.006815 

0.01424 

4 

0.000669 

0.00185 

5 

5.25E-05 

0.00018 

Figure  9  Spatial  Poisson  Process  Probability  Output 

Using  the  spatial  Poisson  process  within  the  model  and  eliminating  the  double  counts  / 
intersections  using  [n*=i(l  —  Vx)  ]*  Pn  *  n  and  then  aggregating  over  all  “n”  the  Risk  Score  = 
.3580 
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Figure  10  contains  the  output  when  using  the  extended  spatial  Poisson  process. 

To 

0  1  2  3  4  5 

OF  0.675366  0.275763  0.04504  0.003678  0.00015  2.45288E-06_ 

1  0  0.730519  0.238627  0.029231  0.001591  3.24894E-05 

From  2  0  0  0.790176  0.193585  0.015809  0.000430334 

3  0  0  0  0.854704  0.139596  0.005699944 

4  0  0  0  0  0.924502  0.075497976 

5) _  0  0  0  0  0  1_ 

Figure  10  Extended  Spatial  Poisson  Process 

Using  the  extended  spatial  Poisson  process  and  subtracting  out  the  double  counts  and  then 


calculating  the  Risk  Score  results  in: 
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=  0.356 

Thus  it  is  clear  that  the  risk  score  using  the  extended  spatial  Poisson  process  is  tighter  because 
the  model  eliminated  any  entity  that  was  removed  when  detennining  the  next  probability. 
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3.7.  Universal  Risk  Score 


In  order  to  standardize  the  risk  score,  the  model  generates  a  scaled  risk  score.  Using  the 
relative  risk  score  from  Table  3  and  removing  deer  and  coyote,  the  relative  risk  score  can  be 
rescaled  using  Equation  12. 

fMJb-a){x-min)+a  (12) 

max-  min 

where  b  =  Rescaled  Max;  a  =  Rescaled  Min;  min  =  Original  Min;  max  =  Original  Max 
This  paper  uses  the  following  factors  for  the  rescale: 
b  =  10  and  a  =  0 

min  =  risk  score  for  entire  AO  if  all  birds  are  Swallows 
max  =  risk  score  for  entire  AO  if  all  birds  are  Vultures 

This  score  is  then  used  to  generate  25/50/25  splits  to  highlight  the  top  25%  in  terms  of  risk  score, 

the  middle  50%  and  the  lowest  25%.  Figure  1 1  illustrates  this  sliding  scale  and  how  a  specific 

scenario  risk  score  relates  to  the  maximum,  75%  and  25%  thresholds. 


MAX  POISSON  WHOLE  AO 
75% 

RiskScore 
25% 

MIN  POISSON  WHOLE  AO 

0  J - 

Figure  1 1  Risk  Score  Scale 

Additionally,  using  a  Risk  Filtering  and  Ranking  Method  (RFRM)  that  compares  the  effect  (total 
risk  score)  and  the  likelihood  of  an  event  happening  (1-P(0)),  the  model  generates  a  graphic  that 
depicts  where  the  scenario’s  risk  level  is  located  on  the  graph  using  a  scale  of 
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low(L)/moderate(M)/high(H)/extremely  high(H+).  This  chart  can  be  tailored  for  different 
components  or  different  airframes,  depending  on  what  is  being  measured. 
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Likelihood  =  (1  -  P(0)) 


Figure  12  Risk  Filtering  and  Ranking  Method 

3.8.  Case  Study 

The  following  scenario  represents  a  hypothetical  flight  through  an  operations  area.  The 
probability  of  hitting  a  bird  is  calculated  and  a  risk  score  is  generated. 

Consider  a  low  level  bombing  range  of  dimension  1km  x  1km  x  3500  ft.  A  B-1B  with 
wings  swept  will  pass  through  the  operations  area.  The  wingspan  of  the  B-1B  is  79ft  and  its 
height  is  26ft.  The  aircraft  will  travel  a  distance  of  1km  through  the  operations  area.  The 
planned  altitude  block  is  1000  -  2000  feet  above  ground  level.  Avian  radar  has  depicted  350 
birds  in  the  operations  area.  Through  historical  data,  specific  bird  species,  their  percent  of  the 
total  and  their  altitude  distribution  have  been  determined  according  to  Table  7. 


38 


Table  7  Scenario  Inputs 


Total  #  birds  in  AO  =  500 

Species 

%  of  Total 

Alpha  "a" 

Beta "P" 

Vultures 

0.1 

8 

50 

Geese 

0.2 

6 

80 

Pelicans 

0.11 

10 

43 

Ducks 

0.14 

5 

151 

Eagles 

0.13 

4.8 

153 

Hawks 

0.09 

4.6 

155 

Gulls 

0.23 

4.2 

159 

The  aircrew  wishes  to  determine  the  probability  of  an  aircraft  strike  and  the  associated 
risk  score  according  to  the  model. 

The  parameters  from  the  above  scenario  are  entered  into  the  model  and  calculations  are 
accomplished  utilizing  A.’ Species Alt,  |A|,  and  |v|  for  each  species.  The  following  figures  (Figure 
13)  highlight  the  calculations  that  are  being  accomplished  in  support  of  the  model  for  vultures 
and  geese. 
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Figure  13  Species  Specific  Calculations 

Given  the  parameters,  the  model  follows  Case  #3  from  Table  6  and  the  bird  strike 
probabilities  and  an  associated  risk  score  can  be  computed.  Figure  14  is  the  final  model  with  all 


the  inputs  and  the  associated  outputs. 
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Calculate  P(0)  *  no  risk  -  Case  #1 


calculate  P(0)  and  risk  full  AO  -  Case  #2 


Calculate  P(O)  and  Risk  (alt  band)  -  Case  #3 
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Figure  14  Aircraft  /  Avian  Encounter  Model 

Note  that  p(0),  the  probability  of  no  bird  strikes,  is  -29%  and  the  probability  of  at  least 
one  bird  strike  (1-  p(0))  is  -71%.  The  risk  score  is  4.54  and  falls  within  the  middle  50%  of  the 
sliding  scale  risk  score  graphic. 

Further  risk  management  fidelity  can  be  gained  from  Figure  15,  the  RFRM  graph.  The 
RFRM  graph  gives  an  overall  factor  of  Moderate  for  this  scenario. 
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Figure  15  RFRM  Graph 

3.9.  Summary 

Using  the  spatial  Poisson  process  to  calculate  a  probability  for  an  aircraft  /  avian 
encounter  provides  a  closed  form  mathematical  model  that  is  fairly  easy  to  compute.  Assessing 
the  risk  is  more  difficult  and  requires  additional  parameter  inputs  and  greater  fidelity  of  those 
parameters.  Additionally,  there  are  several  key  assumptions  that  need  to  be  made  (independence 
of  birds,  small  bird  estimates,  species  percentage,  etc)  when  computing  a  relative  risk  score.  If 
these  assumptions  are  met,  then  planners  can  use  this  model  to  underpin  their  decisions  when 
assessing  overall  risk  management.  The  extended  spatial  Poisson  process  is  a  very  powerful 
method  that  provides  a  higher  degree  of  accuracy  when  determining  encounter  probability  as  it 
allows  for  the  removal  of  entities,  which  are  already  encountered,  from  the  AO. 
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4.  Analysis 


4.1.  Bird  Density 

The  following  scenario  is  for  a  B-1B  transiting  through  an  AO  of  1km  x  1km  x  3500ft.  The 
number  of  total  birds  for  in  the  AO  is  the  independent  variable  in  the  analysis  and  the  total 
number  varies  from  0  to  3000  total  birds.  The  dependent  variable  is  the  probability  of  one  or 
more  bird  strikes  (1-P(0)). 


Figure  16  Sensitivity  Analysis  [1-P(0)] 

As  expected  the  probability  of  one  of  more  strikes  increases  as  the  number  of  total  birds  in  the 
AO  increases.  Figure  16  indicates  that  the  probability  of  one  or  more  bird  strikes  increases  at  a 
lognormal  rate  as  the  number  of  birds  in  the  AO  increases. 
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4.2.  Species  Percentage  /  Risk  Score 


The  following  scenario  is  for  a  B-1B  transiting  through  an  AO  of  1km  x  1km  x  3500ft.  In  this 
analysis,  there  were  7  different  species  present  in  the  AO.  They  consisted  of  vultures,  geese, 
cranes,  pelicans,  ducks,  eagles  and  hawks.  These  7  species  were  identified  by  USAF  as  the  7 
most  damaging  species  for  Air  Combat  Command  (Kelly,  et  al.  1999).  The  independent  variable 
in  this  analysis  was  the  percentage  (of  the  total  #  of  birds)  that  were  vultures.  The  other  six 
species  were  held  at  equal  percentages  of  the  remainder  (calculation:  1/7  *  (1  -  Vulture  %)).  The 
dependent  variable  is  the  Risk  Score.  The  total  number  of  birds  in  the  AO  was  also  varied 
(independently)  from  100  to  500  by  increments  of  100.  The  risk  score  is  scaled  from  0  to  1  for 
brevity. 


Total  Birds  /  %  Vulture  vs  Risk  Score 


0.4  0.6 

Percentage  of  Vultures 


500  birds 
-400  birds 
-300  birds 
200  birds 
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Figure  17  Sensitivity  Analysis  [Risk  Score  /  Vulture  probability] 

Figure  17  highlights  that  as  the  slope  of  the  line  associated  with  the  risk  score  is  steeper  as  the 
total  number  of  birds  in  the  AO  is  increase.  The  risk  score  increases  for  all  level  of  total  number 
of  birds  as  the  percentage  of  vultures  increases  (as  expected). 
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4.3.  Gamma  Altitude  Band 


The  following  scenario  is  for  a  B-1B  transiting  through  an  AO  of  1km  x  1km  x  3500ft  at  an 
altitude  of  300  to  600  meters.  In  this  scenario,  there  are  200  vultures  in  the  AO.  The  independent 
variables  in  this  analysis  are  the  alpha  and  beta  values  associated  with  the  gamma  altitude 
distribution.  They  are  varied  from  .01  to  10  for  the  alpha  value  and  from  1  to  500  for  the  beta 
value.  The  dependent  variable  is  the  risk  score.  The  risk  score  is  scaled  from  0  to  1  for  brevity. 
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From  Figure  18  it  can  be  seen  that  that  the  beta  variable  is  the  most  influential  (50  to  70)  for 
alpha  values  between  4-10.  The  alpha  value  is  the  most  influential  when  alpha  values  are 
below  4  and  they  are  most  influential  from  1.2  -  2.4  for  beta  value  from  100  to  500. 
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5.  Conclusion 


This  paper  combines  the  powerful  methods  of  the  spatial  Poisson  processes  and  the 
extended  spatial  Poisson  processes  with  relative  risk  score  values  and  with  gamma  altitude 
distributions.  The  spatial  Poisson  process  accounts  for  the  spatial  distribution  of  birds  within  a 
bounded  area  of  operations.  The  extended  spatial  Poisson  process  allows  for  over-dispersed  and 
under-dispersed  object  densities  in  this  bounded  space  with  respect  to  the  spatial  Poisson  process 
and  provides  for  a  more  accurate  estimation.  The  relative  risk  scores  account  for  differing 
hazards  associated  with  striking  different  species  of  birds.  The  gamma  altitude  distributions 
provides  further  fidelity  when  determining  the  hazard  associated  with  a  particular  altitude  band 
to  be  flown.  The  model  allows  for  varying  levels  of  user  input  and  provides  numerous  outputs  to 
include  collision  probabilities,  a  raw  risk  score,  a  standardized  and  graphical  relative  risk  score,  a 
gamma  distribution  detennination  tool  and  a  risk  filtering  and  ranking  method  graphic.  It  is 
recommended  that  leaders  and  resident  BASH  /  Wildlife  experts  incorporate  this  model  (or  a 
locally  adjusted  model)  into  their  decision  making  process. 

Future  research  in  the  area  of  aircraft  /  avian  encounters  should  account  for  the  varying 
flock  sizes  (Hemelrijk  and  Hildenbrandt  2011).  The  Shot  noise  Cox  processes  and  the  finite 
Gibbs  point  process  (Moller  and  Waagepeterson  2006)  could  be  utilized  to  model  some  of  the 
clustering  of  birds  within  an  AO.  Additionally,  flock  shapes  and  density  distributions  within  the 
flock  should  be  studied.  This  would  allow  for  a  greater  fidelity  on  types  of  birds  that  are 
observed  on  a  radar  return. 
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7.  Appendix  I  -  Model 


The  complete  model  allows  the  decision  maker  to  input  parameters  for  AO  data,  Aircraft 
data  and  Bird  data.  Any  one  of  the  three  cases  can  be  computed  if  the  parameters  of  species 
percentage  and  gamma  height  distribution  are  known.  Note  that  the  aircraft  data  currently  is 
utilizing  an  ellipse  around  the  aircraft  as  the  encounterable  area.  If  the  decision  maker  knows  the 
exact  frontal  surface  area  or  desires  a  specific  surface  area  (e.g.  engine  intake),  then  the  model 
can  easily  be  modified  to  account  for  this.  The  gamma  distribution  graph  allows  the  decision 
maker  to  match  up  the  specific  bird  height  parameters  with  what  is  known.  The  graph  in  the 
lower  left  corner  is  the  sliding  scale  for  the  risk  score.  The  risk  score  adjusts  based  upon  total 
number  and  species  of  birds  in  the  AO.  The  tables  in  the  lower  right  are  probabilities  P(n)  {n 
equal  >=  1 ,  and  n  =  0,  1,2,34}. 
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The  Risk  Filtering  and  Ranking  Method  is  displayed  via  the  following  graphs  where: 
Y  axis  is  the  effect  (damage)  on  an  aircraft,  increasing 

from :  Marginal  ->  Moderate  ->  S  erious  Critical  ->  Catastrophic : 

X  axis  is  the  likelihood  from  probability  of  a  hit  =  0  to  probability  of  a  hit  =  1 
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The  calculations  for  each  of  the  bird  species  are  computed  utilizing  the  following  tables: 
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Pelican 

Risk  Sum 

16.21572 

Ducks 

Risk  Sum 

11.75228 

14 

Max  Risk 

33.81457 

15 

Max  Risk 

28.43498 

Gamma  Alpha  "a"  = 

10 

Gamma  Alpha  "a"  = 

5 

Gamma  Beta  "(3"  = 

43 

Gamma  Beta  "P"  = 

151 

Gamma  Altitude  Band  % 

0.72143785 

Gamma  Altitude  Band  % 

0.323775914 

Total  of  specific  species  in  AO 

220 

Total  of  specific  species  in  AO 

440 

X' 

5.20381E-07 

X' 

4.67087E-07 

Risk  Multiplier 

44 

Risk  Multiplier 

37 

MAX 

MAX 

n 

P(n) 

Risk  P(n) 

Risk 

n 

P(n) 

Risk 

P(n)  Risk 

0 

0.431933625 

0  0.22 

0 

0 

0.470712534 

0 

0.22  0 

1 

0.362601086 

9.063185  0.33 

11.43 

1 

0.35468552 

6.946032 

0.33  9.61 

2 

0.152198787 

4.849582  0.25 

11.53671 

2 

0.133629136 

3.377504 

0.25  9.701325 
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0.042589449 

1.72576  0.13 

6.528561 
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0.033563528 

1.102446 

0.13  5.489926 
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0.008938283 

0.462348  0.05 

2.874446 
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0.267607 

0.05  2.417147 
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1.035169 
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0.01  0.870483 
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0.308747 
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0.00  0.259628 
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7.65367E-09 

8.04E-07 

0.00  0.000448 
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1.57904E-09 
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8.07E-05 

11 

5.24282E-10 

6.06E-08 

0.00  6.78E-05 
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1.10465E-10 

1.7E-08  0.00 

1.11E-05 

12 

3.29209E-11 

4.15E-09 

0.00  9.33E-06 
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7.13335E-12 

1.19E-09  0.00 

1.4E-06 
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2.61E-10 
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7.66E-11  0.00 

1.63E-07 
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1.76E-08 

15 

5.15907E-15 
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0.00  1.26E-11 

19 
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3.  HE- 17  0.00 
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3.57E-18 

0.00  1.06E-12 
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5.36456E-21 

1.37E-18  0.00 

IE- 13 

20 

6.73574E-22 

1.41E-19 

0.00  8.42E-14 
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Risk  Sum 
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Risk  Sum 
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16 

Max  Risk 

23.8239 
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Max  Risk 

19.21282 

Gamma  Alpha  "a"  = 

4.8 

Gamma  Alpha  "a"  = 

4.6 

Gamma  Beta  "P"  = 
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Gamma  Beta  "P"  = 
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Gamma  Altitude  Band  % 

0.340323459 

Gamma  Altitude  Band  % 

0.356522721 

Total  of  specific  species  in  AO 
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Total  of  specific  species  in  AO 
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X' 

8.92652E-07 

X' 
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Risk  Multiplier 
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Risk  Multiplier 
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MAX 
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P(n) 

Risk 
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Risk 
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Risk 
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Risk 
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1 
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1 
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2 
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4 
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15 
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IE-08 

16 
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16 
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17 
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19 
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19 
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20 
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20 
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Risk  Sum 
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Risk  Sum 

2.292915 

18 

Max  Risk 

18.44431 

19 

Max  Risk 

16.90728 

Gamma  Alpha  "a"  = 

4.4 

Gamma  Alpha  "a"  = 

4.2 

Gamma  Beta  "(3"  = 

141 

Gamma  Beta  "P"  = 

159 

Gamma  Altitude  Band  % 

0.427237501 

Gamma  Altitude  Band  % 

0.386576789 

Total  of  specific  species  in  AO 

160 

Total  of  specific  species  in  AO 

X' 

2.24125E-07 

X' 

2.28144E-07 

Risk  Multiplier 

24 

Risk  Multiplier 

22 

MAX 

MAX 

n 

P(n) 

Risk  P(n) 

Risk 

n 

P(n) 

Risk 

P(n)  Risk 

0 

0.696589117 

0  0.22 

0 

0 

0.692087305 

0 

0.22  0 

1 

0.251858444 

1.833998  0.33 

6.24 

1 

0.254718004 

1.72548 

0.33  5.72 

2 

0.045530912 

0.496092  0.25 

6.292751 

2 

0.046873611 

0.473292 

0.25  5.768355 

3 

0.005487379 

0.0856  0.13 

3.561033 

3 

0.005750504 

0.083013 

0.13  3.26428 

4 

0.000496004 

0.01026  0.05 

1.567879 

4 

0.000529108 

0.010126 

0.05  1.437223 

5 

3.5867E-05 

0.000927  0.01 

0.564638 

5 

3.89469E-05 

0.000931 

0.01  0.517584 

6 

2.16134E-06 

6.7E-05  0.00 

0.168408 

6 

2.38903E-06 

6.85E-05 

0.00  0.154374 

7 

1.11636E-07 

4.04E-06  0.00 

0.04232 

7 

1.25609E-07 

4.2E-06 

0.00  0.038793 

8 

5.04539E-09 

2.09E-07  0.00 

0.009142 

8 

5.7787E-09 

2.21E-07 

0.00  0.00838 

9 

2.0269E-10 

9.43E-09  0.00 

0.001729 

9 

2.36313E-10 

1.02E-08 

0.00  0.001585 

10 

7.32845E-12 

3.79E-10  0.00 

0.000291 

10 

8.69732E-12 

4.16E-10 

0.00  0.000267 

11 

2.40879E-13 

1.37E-11  0.00 

4.4E-05 

11 

2.90999E-13 

1.53E-11 

0.00  4.03E-05 

12 

7.25768E-15 

4.5E-13  0.00 

6.05E-06 

12 

8. 92502 E- 15 

5.12E-13 

0.00  5.55E-06 

13 

2.01853E-16 

1.36E-14  0.00 

7.63E-07 

13 

2.52676E-16 

1.57E-14 

0.00  7E-07 

14 

5.21298E-18 

3.77E-16  0.00 

8.89E-08 

14 

6.64256E-18 

4.45E-16 

0.00  8.15E-08 

15 

1. 25653 E- 19 

9.74E-18  0.00 

9.61E-09 

15 

1.62983E-19 

1.17E-17 

0.00  8.8E-09 

16 

2.83945E-21 

2.35E-19  0.00 

9.69E-10 

16 

3.74905E-21 

2.87E-19 

0.00  8.88E-10 

17 

6.039E-23 

5.31E-21  0.00 

9.17E-11 

17 

8.11655E-23 

6.6E-21 

0.00  8.4E-11 

18 

1.21303E-24 

1.13E-22  0.00 

8.16E-12 

18 

1.65958E-24 

1.43E-22 

0.00  7.48E-12 

19 

2.30834E-26 

2.27E-24  0.00 

6.86E-13 

19 

3.21472E-26 

2.92E-24 

0.00  6.29E-13 

20 

4. 173E-28 

4.31E-26  0.00 

5.46E-14 

20 

5.91577E-28 

5.66E-26 

0.00  5.01E-14 

Herons 

Risk  Sum 

5.028598 

Mourning  Dove 

Risk  Sum 

6.170355 

20 

Max  Risk 

16.90728 

21 

Max  Risk 

13.06472 

Gamma  Alpha  "a"  = 

4 

Gamma  Alpha  "a"  = 

3.8 

Gamma  Beta  "P"  = 

141 

Gamma  Beta  "P"  = 

141 

Gamma  Altitude  Band  % 

0.454207518 

Gamma  Altitude  Band  % 

0.461588324 

Total  of  specific  species  in  AO 

250 

Total  of  specific  species  in  AO 

340 

X' 

3.72301E-07 

X' 

5.14557E-07 

Risk  Multiplier 

22 

Risk  Multiplier 

17 

MAX 

MAX 

n 

P(n) 

Risk 

P(n) 

Risk 

n 

P(n) 

Risk 

P(n) 

Risk 

0 

0.548482842 

0 

0.22 

0 

0 

0.436010851 

0 

0.22 

0 

1 

0.329418401 

3.272237 

0.33 

5.72 

1 

0.36192744 

3.470094 

0.33 

4.42 

2 

0.098924227 

1.317896 

0.25 

5.768355 

2 

0.150215839 

1.837957 

0.25 

4.457366 

3 

0.019804607 

0.356613 

0.13 

3.26428 

3 

0.041564129 

0.648244 

0.13 

2.522398 

4 

0.002973658 

0.06998 

0.05 

1.437223 

4 

0.008625473 

0.171911 

0.05 

1.110581 

5 

0.000357195 

0.010476 

0.01 

0.517584 

5 

0.001431981 

0.035368 

0.01 

0.399952 

6 

3.57552E-05 

0.001258 

0.00 

0.154374 

6 

0.000198112 

0.005863 

0.00 

0.119289 

7 

3.06779E-06 

0.000126 

0.00 

0.038793 

7 

2.34929E-05 

0.000811 

0.00 

0.029977 

8 

2.30314E-07 

1.08E-05 

0.00 

0.00838 

8 

2.43765E-06 

9.62E-05 

0.00 

0.006475 

9 

1.53696E-08 

8.11E-07 

0.00 

0.001585 

9 

2.24829E-07 

9.98E-06 

0.00 

0.001225 

10 

9.23099E-10 

5.41E-08 

0.00 

0.000267 

10 

1.86628E-08 

9.2E-07 

0.00 

0.000206 

11 

5.04011E-11 

3.25E-09 

0.00 

4.03E-05 

11 

1.40834E-09 

7.64E-08 

0.00 

3.12E-05 

12 

2.52257E-12 

1.77E-10 

0.00 

5.55E-06 

12 

9.74206E-11 

5.77E-09 

0.00 

4.29E-06 

13 

1. 16543 E- 13 

8.88E-12 

0.00 

7E-07 

13 

6.22059E-12 

3.99E-10 

0.00 

5.41E-07 

14 

4.99968E-15 

4.  IE- 13 

0.00 

8.15E-08 

14 

3.68831E-13 

2.55E-11 

0.00 

6.29E-08 

15 

2.00187E-16 

1.76E-14 

0.00 

8.8E-09 

15 

2.04108E-14 

1.51E-12 

0.00 

6.8E-09 

16 

7.51451E-18 

7.05E-16 

0.00 

8.88E-10 

16 

1.05892E-15 

8.36E-14 

0.00 

6.86E-10 

17 

2.65483E-19 

2.65E-17 

0.00 

8.4E-11 

17 

5.17059E-17 

4.33E-15 

0.00 

6.49E-11 

18 

8.85826E-21 

9.35E-19 

0.00 

7.48E-12 

18 

2.38447E-18 

2.12E-16 

0.00 

5.78E-12 

19 

2.80014E-22 

3.12E-20 

0.00 

6.29E-13 

19 

1.04175E-19 

9.76E-18 

0.00 

4.86E- 13 

20 

8.40881E-24 

9.86E-22 

0.00 

5.01E-14 

20 

4.32371E-21 

4.26E-19 

0.00 

3.87E-14 

Owls 

Risk  Sum 

8.654626 

American  Kestrel 

Risk  Sum 

8.533788 

22 

Max  Risk 

12.29621 

23 

Max  Risk 

10.75918 

Gamma  Alpha  "a"  = 

3.6 

Gamma  Alpha  "a"  = 

3.4 

Gamma  Beta  "(3"  = 

141 

Gamma  Beta  "P"  = 

141 

Gamma  Altitude  Band  % 

0.46417771 

Gamma  Altitude  Band  % 

0.461544896 

Total  of  specific  species  in  AO 

460 

Total  of  specific  species  in  AO 

510 

X' 

7.00071E-07 

X' 

7.71764E-07 

Risk  Multiplier 

16 

Risk  Multiplier 

14 

MAX 

MAX 

n 

P(n) 

Risk 

P(n) 

Risk 

n 

P(n) 

Risk 

P(n) 

Risk 

0 

0.323239902 

0 

0.22 

0 

0 

0.287936567 

0 

0.22 

0 

1 

0.365054378 

3.952868 

0.33 

4.16 

1 

0.358485367 

3.5737 

0.33 

3.64 

2 

0.206138997 

2.834532 

0.25 

4.195168 

2 

0.223159843 

2.854298 

0.25 

3.670772 

3 

0.077601747 

1.270657 

0.13 

2.374022 

3 

0.092612456 

1.380307 

0.13 

2.077269 

4 

0.021910087 

0.441223 

0.05 

1.045253 

4 

0.028825976 

0.519783 

0.05 

0.914596 

5 

0.004948877 

0.121846 

0.01 

0.376425 

5 

0.007177755 

0.157121 

0.01 

0.329372 

6 

0.000931511 

0.027385 

0.00 

0.112272 

6 

0.001489402 

0.038843 

0.00 

0.098238 

7 

0.000150287 

0.00515 

0.00 

0.028213 

7 

0.000264904 

0.008048 

0.00 

0.024687 

8 

2.12161E-05 

0.000831 

0.00 

0.006095 

8 

4.12262E-05 

0.001431 

0.00 

0.005333 

9 

2.66229E-06 

0.000117 

0.00 

0.001153 

9 

5.70303E-06 

0.000223 

0.00 

0.001009 

10 

3.00668E-07 

1.47E-05 

0.00 

0.000194 

10 

7.10035E-07 

3.08E-05 

0.00 

0.00017 

11 

3.08694E-08 

1.66E-06 

0.00 

2.93E-05 

11 

8.03641E-08 

3.84E-06 

0.00 

2.57E-05 

12 

2.90522E-09 

1.71E-07 

0.00 

4.04E-06 

12 

8.33787E-09 

4.34E-07 

0.00 

3.53E-06 

13 

2.52388E-10 

1.61E-08 

0.00 

5.09E-07 

13 

7.98521E-10 

4.5E-08 

0.00 

4.45E-07 

14 

2.03598E-11 

1.4E-09 

0.00 

5.92E-08 

14 

7.10122E-11 

4.31E-09 

0.00 

5.18E-08 

15 

1.5329E-12 

1.13E-10 

0.00 

6.4E-09 

15 

5. 89408 E- 12 

3.84E-10 

0.00 

5.6E-09 

16 

1.082E-13 

8.47E-12 

0.00 

6.46E-10 

16 

4.58639E-13 

3.18E-11 

0.00 

5.65E-10 

17 

7.18804E-15 

5.98E-13 

0.00 

6.11E-11 

17 

3.3589E-14 

2.48E-12 

0.00 

5.35E-11 

18 

4.50994E-16 

3.97E-14 

0.00 

5.44E-12 

18 

2.32327E-15 

1.81E-13 

0.00 

4.76E-12 

19 

2.68071E-17 

2.49E-15 

0.00 

4.57E-13 

19 

1.52237E-16 

1.25E-14 

0.00 

4E-13 

20 

1.51374E-18 

1.48E-16 

0.00 

3.64E-14 

20 

9.47686E-18 

8.22E-16 

0.00 

3.19E-14 
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Sparrows 


Blackbirds  -  Starling 

Risk  Sum 

9.274724 

26 

Max  Risk 

6.916616 

Gamma  Alpha  "a"  = 

2.8 

Gamma  Beta  "3"  = 

141 

Gamma  Altitude  Band  % 

0.420114193 

Total  of  specific  species  in  AO 

880 

X' 

1.21213E-06 

Risk  Multiplier 

9 

MAX 

n 

P(n) 

Risk 

P(n) 

Risk 

0 

0.141504743 

0 

0.22 

0 

1 

0.276701493 

2.137922 

0.33 

2.34 

2 

0.2705341 

3.023779 

0.25 

2.359782 

3 

0.176336114 

2.15658 

0.13 

1.335387 

4 

0.086202881 

1.157803 

0.05 

0.587955 

5 

0.033712603 

0.517208 

0.01 

0.211739 

6 

0.010987061 

0.195453 

0.00 

0.063153 

7 

0.003069192 

0.062999 

0.00 

0.01587 

8 

0.000750196 

0.017544 

0.00 

0.003428 

9 

0.000162994 

0.004285 

0.00 

0.000648 

10 

3. 18723E-05 

0.000931 

0.00 

0.000109 

11 

5.66579E-06 

0.000182 

0.00 

1.65E-05 

12 

9.23252E-07 

3.24E-05 

0.00 

2.27E-06 

13 

1.38873E-07 

5.27E-06 

0.00 

2.86E-07 

14 

1.93968E-08 

7.93E-07 

0.00 

3.33E-08 

15 

2.52859E-09 

1.11E-07 

0.00 

3.6E-09 

16 

3.09029E-10 

1.44E-08 

0.00 

3.63E-10 

17 

3.5546E-11 

1.76E-09 

0.00 

3.44E-11 

18 

3.86153E-12 

2.03E-10 

0.00 

3.06E-12 

19 

3.97417E-13 

2.21E-11 

0.00 

2.57E-13 

20 

3.88559E-14 

2.27E-12 

0.00 

2.05E-14 

Risk  Sum  0.392892 


27 

Max  Risk 

3.074052 

Gamma  Alpha  "a"  = 

2.6 

Gamma  Beta  "P"  = 

141 

Gamma  Altitude  Band  % 

0.395301966 

Total  of  specific  species  in  AO 

170 

X' 

2.20332E-07 

Risk  Multiplier 

4 

MAX 

n 

P(n) 

Risk 

P(n) 

Risk 

0 

0.700863803 

0 

0.22 

0 

1 

0.249116222 

0.298079 

0.33 

1.04 

2 

0.044273147 

0.079556 

0.25 

1.048792 

3 

0.005245507 

0.013513 

0.13 

0.593506 

4 

0.000466118 

0.001593 

0.05 

0.261313 

5 

3.31356E-05 

0.000141 

0.01 

0.094106 

6 

1.96296E-06 

1.01E-05 

0.00 

0.028068 

7 

9.9674E-08 

5.96E-07 

0.00 

0.007053 

8 

4.42854E-09 

3.02E-08 

0.00 

0.001524 

9 

1.74898E-10 

1.34E-09 

0.00 

0.000288 

10 

6.21662E-12 

5.31E-11 

0.00 

4.85E-05 

11 

2.00877E-13 

1.89E-12 

0.00 

7.33E-06 

12 

5.95E-15 

6.1E-14 

0.00 

1.01E-06 

13 

1.62683E-16 

1.81E-15 

0.00 

1.27E-07 

14 

4. 13031E-18 

4.94E-17 

0.00 

1.48E-08 

15 

9.78723E-20 

1.25E-18 

0.00 

1.6E-09 

16 

2.17424E-21 

2.97E-20 

0.00 

1.62E-10 

17 

4.54598E-23 

6.6E-22 

0.00 

1.53E-11 

18 

8.97684E-25 

1.38E-23 

0.00 

1.36E- 12 

19 

1.67934E-26 

2.72E-25 

0.00 

1.14E- 13 

20 

2.98453E-28 

5.1E-27 

0.00 

9.  HE- 15 

Swallows  Risk  Sum 

28  Max  Risk 
Gamma  Alpha  "a"  = 

Gamma  Beta  "(3"  = 

Gamma  Altitude  Band  %  0.365606619 
Total  of  specific  species  in  AO^|  880 

X'  1.05486E-06 
Risk  Multiplier  2 


1.760176 

1.537026 


n  P(n) 

Risk 

P(n) 

Risk 

0  0.182370222 

0 

0.22 

0 

1  0.310342411 

0.50749 

0.33 

0.52 

2  0.264057396 

0.595592 

0.25 

0.524396 

3  0.149783607 

0.372949 

0.13 

0.296753 

4  0.063722308 

0.179864 

0.05 

0.130657 

5  0.02168746 

0.071644 

0.01 

0.047053 

6  0.006150985 

0.023855 

0.00 

0.014034 

7  0.001495319 

0.006724 

0.00 

0.003527 

8  0.000318076 

0.001632 

0.00 

0.000762 

9  6.01417E-05 

0.000347 

0.00 

0.000144 

10  1.02344E-05 

6.56E-05 

0.00 

2.42E-05 

11  1.58328E-06 

1.12E-05 

0.00 

3.67E-06 

12  2.24524E-07 

1.73E-06 

0.00 

5.04E-07 

13  2.93905E-08 

2.45E-07 

0.00 

6.36E-08 

14  3.57245E-09 

3.21E-08 

0.00 

7.4E-09 

15  4.05287E-10 

3.9E-09 

0.00 

8E-10 

16  4.31052E-11 

4.42E-10 

0.00 

8.08E-11 

17  4.31487E-12 

4.7E-11 

0.00 

7.64E-12 

18  4.07927E-13 

4.71E-12 

0.00 

6.8E-13 

19  3.65356E-14 

4.45E-13 

0.00 

5.72E-14 

20  3.10866E-15 

3.99E-14 

0.00 

4.55E-15 
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8.  Appendix  II  -  Numerical  Matrix  Exponentiation  Code 


The  following  is  VBA  code  for  numerical  calculation  of  the  Extended  Spatial  Poisson  Process: 

Opt  ion  Explicit  fj 

Public  Sub  EXPO() 

Dim  iRow  As  Integer,  iColumn  As  Integer,  ilteration  As  Integer 

Dim  i  As  Long,  N  As  Long 

Dim  V  As  Long 

Dim  Lambda  As  Double 

Dim  vDistance  As  Double 

Dim  vLength  As  Double 

Dim  vHeight  As  Double 

Dim  aircraftA  As  Double 

Dim  arrayA()  As  Variant  'A  matrix 

Dim  arrayl()  As  Variant  'I  Matrix 

Dim  arrayQO  As  Variant 

Dim  arrayPower ()  As  Variant 

Dim  arrayPowerFact ( )  As  Variant 

Application . ScreenUpdating  =  False 

' Input 

N  =  Sheetl. Range ("Bl”)  '#  of  birds  in  u(S) 

V  =  Sheetl. Range ("B2")  'volume  of  AO 

ilteration  =  Sheetl . Range ( "B3")  'Iteration  for  series 
vDistance  =  Sheetl . Range ( "B4" )  'distance  aircraft  travels  in  AO 
vLength  =  Sheetl . Range ( "B5" )  'wingspan 
vHeight  =  Sheetl . Range ( "B 6" )  'height 

'Initialize  Variables 
Lambda  =  -1  *  Log(l  -  (1  /  V) ) 

aircraftA  =0.5  *  vLength  *  0.5  *  vHeight  *  vDistance  *  3.14 

Sheet2 . Select 
Clear_Sheet 

'Build  Identity  Matrix 

ReDim  arrayl(l  To  N  +  1,  1  To  N  +  1) 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

If  iRow  =  iColumn  Then 

arrayl (iRow,  iColumn)  =  1 

Else 

arrayl (iRow,  iColumn)  =  0 
End  If 

Next 

Next 


'Build  A  Matrix 


ReDim  arrayA(l  To  N  +  1,  1  To  N  +  1) 
For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

If  iRow  =  iColumn  Then 

arrayA(iRow,  iColumn)  =  -1 
Elself  iColumn  =  iRow  +  1  Then 
arrayA(iRow,  iColumn)  =  (N 
Elself  iColumn  =  iRow  +  1  Then 
Else 

arrayA(iRow,  iColumn)  =  0 
End  If 

Next 

Next 


(N  -  (iRow  -  1))  *  Lambda  *  aircraftA 
(iRow  -  1) )  *  Lambda  *  aircraftA 


' Initial  summation  I  +  A 

ReDim  arrayQ(l  To  N  +  1,  1  To  N  +  1) 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

arrayQ(iRow,  iColumn)  =  arrayl (iRow,  iColumn)  +  arrayA(iRow,  iColumn) 

Next 

Next 
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'Initial  step  matrix  calculation  —  A  squared 
ReDim  arrayPower (1  To  N  +  1,  1  To  N  +  1) 

arrayPower ()  =  Application. WorksheetFunction .MMult (arrayA,  arrayA) 


' divide  by  2 ! 

ReDim  arrayPowerFact ( 1  To  N  +  1,  1  To  N  +  1) 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

arrayPowerFact (iRow,  iColumn)  =  arrayPower (iRow,  iColumn)  /  Application . WorksheetFunction . Fact 

(2) 

Next 

Next 

’Add  I  +  A  +  (A*2) / 2 ! 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

arrayQ(iRow,  iColumn)  =  arrayQ(iRow,  iColumn)  +  arrayPowerFact (iRow,  iColumn) 

Next 

Next 

' Subsequent  step  matrix  calculations 
For  i  =  3  To  ilteration 

'take  next  power  of  A 

arrayPower ()  =  Application. WorksheetFunction .MMult (arrayPower ,  arrayA) 

' divide  by  i ! 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

arrayPowerFact (iRow,  iColumn)  =  arrayPower (iRow,  iColumn)  /  Application. WorksheetFunction. 

Fact  (i) 

Next 

Next 

' add  to  Q 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

arrayQ(iRow,  iColumn)  =  arrayQ(iRow,  iColumn)  +  arrayPowerFact (iRow,  iColumn) 

Next 

Next 


Next 

For  iRow  =  1  To  N  +  1 

For  iColumn  =  1  To  N  +  1 

Sheet2 . Cells (iRow,  iColumn)  =  arrayQ(iRow,  iColumn) 

Next 

Next 

Application . ScreenUpdating  =  True 
Sheet2 . Select 

End  Sub 


58 


9.  Appendix  III  -  Road  Show 


Sjf  Air  Force  Institute  of  Technology  MU 


1  The  AFITof  Today  is  die  Air  Force  of  Tomorrow* 


ESTIMATING  BIRD/ AIRCRAFT 
COLLISION  PROBABILITIES  AND 
RISK  UTILIZING  SPATIAL 
POISSON  PROCESSES 

A  Graduate  Research  Project  by: 

Maj  Brady  Vaira 
Dr.  Jeffery  Cochran,  Advisor 

Sponsor:  AMC  SE/SEF 

30  May  2012 

ah  university.  The  Intellectual  and  Leadership  Center  or  the  Air  Force 
Aim  High.  .Fly- Fight- Win 

W  Overview  <-AFi[ 

The  AFITof  Today  is  the  Air  Force  of  Tom orrow^^^^^^^^^ 


•  Background 

•  Problem  Statement 

•  Research  Objectives 

•  Literature  Review 

•  Methodology 

•  Analysis 

•  Conclusion 

•  Future  Research 

•  References 

— i  Air  university.  The  Intellectual  and  Leadership  Center  or  the  Air  Force 
Aim  High. ..Fly- Fight-  Win 
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V  Background  f^AFlJ 

The  AFIT  of  Today  is  the  Air  Force  of  Tomon-ow^^^^^^^^^ 

Random  collisions  of  relatively 
small  numbers  of  entities  in  very 
large,  yet  bounded,  spaces  are 
rare,  but  not  impossible! 


15 Jan2009  3  Feb  2009 


10  Feb  2009 


US  Airways  flight  1549  HMSVanguard&  Iridium  33  and  Kosmos  2251 

LeTriomphant 

Air  unn^rsity.  The  intellectual  and  Leadership  Center  of  the  Air  Force 

Aim  High. Fly  -  Fight-  Win  3 


W  Background  fjFAFIT 

TheAFITof  Today  is  the  Air  Force  of  Tomorrow^^^^^^^^^ 

□  Federal  Aviation  Administration  (FAA)data  from1990  -  2010 

□  121,000  (civil  and  U  S  Air  Force)  wildlife  strikes 

□  United  States  Air  Force  (USAF)  data  from  1 985  -  201 1 

□  95,383  wildlife  strikes 

□  33  fatalities 

□  39  aircraft  lost 

□  Damages  >  $820M 

□  Stakeholders 

□  FAA 

□  DoD 

□  NASA 

□  DOE  (specifically  wind  energy) 

□  Dept  of  Fish  and  Game/Wildlife 

*  U  S  Air  Force  Safety  Center  (2012) 

Air  university:  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
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W  Background  4^AFjT 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow 

□  Federal  Aviation  Administration  (FAA) 

□  92%  of  the  bird  strikes  to  commercial  aircraft  occur  at 
or  below  3,500  ft  above  ground  level  (AGL) 


□  United  States  Air  Force  (USAF) 

□  96.72%  of  its  bird  strikes  below  3,500  ft  AGL. 


USAF  Wildife  Strike; 
Phase  of  Flight 

by  Phase  of  Operation  (1995  -  2011) 

Cost  %  of  Total  Count 

4*  of  Total 

Takeoff  Initial  Climb 

S13T.035  238 

32.131* 

7299 

12254* 

Enroute  Air  Work  Air-to-Air  Air  Refueling 

S15.1S3S8C 

4.504* 

2830 

4.754* 

Flight  Demonstration 

Sl,8'8,791 

0.444* 

27 

0.054* 

Lo  tv  Level  Air-to-Ground  Air  Delivery 

S1T4.151930 

40.8-4* 

<5949 

11.664* 

Hover 

SO 

0.004* 

10 

0.024* 

TVaffc  Pattern  Go-Around 

S31 .446.708 

7384* 

7618 

12.784* 

Initial  Approach  Final  Approach  Landin' 

S46. 825352 

10.994* 

16048 

26924* 

Parked  Ground  Ops 

S3 .561. 521 

0.844* 

423 

0.714* 

Unlawwn 

S12.079.198 

2.834* 

18400 

30.874* 

*U  S  Atr  Force  Safety  Center  (2012) 

Air  university.  The  intellectual  and  Leadership  Center  of  the  Air  Force 

Aim  High.  .Fly- Fight- win 5 


Background 


The  AFIT  of  Today  is  die  Air  Force  ofTomorrowt^^ 

□  Current  Advisory  Systems 

Avian  Hazard  Advisory  System  (AHAS) 


U1AS  KINK  FOR  VRMJ 

VMC4JX  •  JOW/Ot/ll  MW 
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HrXBAO 

MX 

WMMC  Mll/M/II  MW 
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HI  XKAD 
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IOW 

IOW 

hioao 

MX 

«04«  r  NW/os/imu 
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IOW 

M1IAO 

MX 

MOOf  U  MD/M/tl  } 

NOMMtr 

HI  XIAO 

MX 

H  l«  MW 

IOW 

IOW 

MIXEAO 

MX 

mum  i  mi/«/ la 11U 

MVtii 

MWHf 

MIRAO 

MX 

vmoxji  i  fou/oi/ia  tut 

wvm 

WVfU 

HrXBAO 

MX 

««00>  K  M12/M/1I  liUfNOMMt! 

MOOtXAtt 

hi  brad 

MX 

VliMJt  t  NO/OS/IIMV  HOOIMir 

MOWIATI 

HIXRAO 

MX 

HHX  n  Mii/n/UMit 

IOW 

IOW 

HIXNAO 

MX 

WMJH  N  tOli/m/  IS  MW 

tow 

tow 

HI  XIAO 

MX 

WMW  O  M1WM/II  MW 

IOW 

IOW  Xf  IKAO 

MX 

vk«w.. R«.«  V:.\T  *.\1'  - 

H  WARDS 

DAMS: 

VROtiE-F  BEAVER  CREEK  DAM 

VRMJE  F  CHERRYSTONE  CREEK  DAM  -1 

*  United  States  AirForce  BASH  (2012) 

Air  university .  The  intellectual  and  Leadership  Center  or  me  AirForce 
Ann  High. ..Fly- Fight  ■  Win 
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W  Background  -  .-.-U 1 1 

The  AFIT  of  Tod  ayis  the  Air  Force  ofTomorrow^^^^^^^^^- 

□  CurrentAdvisory  Systems 

Bird  Avoidance  Model  (BAM) 


*  United  States  Air  Force  RASH  (2012) 

i  Air  university:  The  intellectual  &id  Leadership  Center  of  me  Air  Force 

Aim  High-Fly  -  Fight  -  win  7 


Problem  Statement 

■  The  AFITof  Today  is  die  Air  Force  of  Tomorrows 


•  The  purpose  of  this  research  is  to  create  a  closed  form 
mathematical  model  to  determine  the  probability  of  an 
aircraft/wildlife  collision  encounter  (strike)  in  a  defined 
space  (e  g  on  a  segment  of  a  Low  Level  route/Air-to- 
Ground/Air  Delivery)  with  a  set  of  given  parameters 


•  Generate  a  risk  score  associated  with  the  encounter(s). 

•  Develop  a  Decision  Support  System  (DSS)  in  order  to 
underpin  decisions  by  pilots  and  planners  when  coupled 
with  existing  alert  and  advisory  systems 

— — Aw  University  Tne  Intellectual  and  Leadership  Center  or  the  A  w  Force 

Aim  High. ..Fly  •  Fight  ■  Win  8 
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Research  Objectives 

The  AFITof  Today  is  die  Air  Force  of  Tomorrow"^ 


•  Study  the  systems  that  provide  avian  detection,  alerts  and 
advisories;  evaluate  the  level  of  fidelity  they  provide. 


•  Determine  avian  characteristics,  geographical  and 
environmental  conditions,  and  other  factors  that  drive  avian 
population  densities  and  spatial  patterns. 

•  Conduct  research  on  Spatial  Poisson  Processes  and  relate 
them  to  aircraft/avian  encounters  for  a  defined  space. 

•  Highlight  and  compare  existing  aircraft/avian  collision 
mathematical  models. 

•  Conduct  research  on  Extended  Spatial  Poisson  Processes 
and  how  they  relate  to  encounter  models. 

•  I  mplement  techniques  for  performing  matrix  exponentiation. 


Air  university:  The  Intellectual  and  Leadership  Center  of  the  Air  Force 

Aim  High.  .Fly- Fight- Wm  9 


Literature  Review 

The  AFITof  Today  is  die  Air  Force  of  Tomorrow'. 


•  Avian  Characteristics  and  Environmental  factors  that 
influence  populations 

•  Land  use  practices 

•  Water  management  facilities 

•  Agricultural  activities 


•  Avian  movements  can  be  partitioned  into  three 
categories;  migrating  birds,  commuting  birds  and 
resident  birds. 

•  During  peak  migration  in  May  and  September  densities 
can  be  as  high  as  1 500  or  more  birds  per  cubic 
kilometer. 

•  FAA  (2010) 

Air  University  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
Atm  High. ..Fly -Fight  -  Win 
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Literature  Review 


The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow' 


Land  use  practices  -fr 

•  Waste  disposal  operations 

•  Trash  Transfer  Stations 
■  Recycling  centers 

•  Construction 

•  GolfCourses 

Water  managementf  acilities^ 

•  Storm  water  management 
facilities 

•  Wastewatertreatment 
facilities 

•  Artificial  marshes 

•  Wastewater  discharge  and 
sludge  disposal 

Agricultural  activities ☆ 

•  Crop  protection 

•  Livestock  production 

•  Aquaculture 


Aw  university  The  intellectual  and  Leadership  Center  of  the  Air  Force 

Aim  High.  .Fly- Fight- Wm  11 


Literature  Review 

The  AFITof  Today  is  the  Air  Force  of  Tomorrow'. 


•  Avian  Altitude  Density  Distributions 

•  Various  literatures  all  reported  that  the  gamma  distribution 
is  good  fit  for  bird  height  distribution. 


*  Shamoun-Baranes,  et  al  (2005)  &  Stumpf  et  al  (2011) 

— Air  University.  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
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'w*  Literature  Review  <-AHl 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow^^^^^^^^*- 

•  Avian  Detection  and  Advisory 

•  Radar  has  been  used  extensively  to  monitor  bird 
movements  and  warn  the  relevant  personnel 

•  Radar  used  to  monitor  bird  movements  in  the  U  S.  is  the 


WSR-88DN EXRAD  radar. 


Sg?  Literature  Review  -VAN  1 

The  AF1T  of  Today  is  the  Air  Force  of  Tomorrow^^^^^^^^^ 


•  Risk  Level  Ranking 

•  Criteria 

•  Damage 

•  MajorDamage 

•  Effect  on  Flight 

•  All  pairwise  ranked 

•  Composite  score 

•  Strongly  correlated  to 
mass  of  bird 
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*  Colbeer  Wnght  and  Cleary  (2000) 
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Literature  Review 


^AFIT 


1  The  AFITof  Today  is  die  Air  Force  of  Tomorrow? 


Spatial  Poisson  Process 

•  One  of  the  simplest  and 
fundamental  spatial  point 
processes  is  the  completely 
random,  or  Spatial  Poisson  Point 
Process. 

•  It  possesses  the  property  of  '‘no 
interaction"  between  points  or 
“complete  spatial  randomness" 

p„  -  .-WI .  PL-  '*!>• 

Fn  n! 


*  t.falt&r  and  Waagepet  arson  (2006) 

Air  university  .  The  intellectual  ana  Leadership  CenteroT  the  Air  Force 
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Literature  Review 

The  AFITof  Today  is  die  Air  Force  of  Tomorrow'. 


•  Extended  Spatial  Poisson  Process 

•  Based  upon  generalizing  the  simplest  Markov  birth 
process,  the  Poisson  process. 


•  ESPP's  involve  representing  a  discrete  distribution  as 
the  distribution  of  the  number  of  events  occurring  in  a 
finite  time  interval  of  a  state-dependent  Markov  birth- 
death  process. 

•  Holzmannand  Cochran  (2012)  furthered  Faddy's 
studies  of  spatial  data  with  respect  to  Extended  SPP 

•  They  utilize  a  linearly  decreasing  arrival  rate  for  the 
birth-death  process  of  the  Markov  transition  in  order  to 
remove  encounter-able  elements  from  the  space. 

*  Podhch  Faddy  and  Smyth  (1999)  &  Holzmann  and  Cochran  (2012) 

Air  University.  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
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Literature  Review 

1  TheAFITofTodayis  the  Air  Force  of  Tomorrow* 


Extended  Spatial  Poisson  Process 


Q  = 


[-NX  NX  0 

0  -(.N  -  1)A  (N  -  1)A 

0  0  -(N  -  2)A 


0 


0 


...  On 

0 

0 

0  OJ 


Where  A  =  In  (1  -  )  and  X  =  #  of  entities  in  bounded  region 

This  then  leads  to  die  probability  distribution  for  the  number  of  encounters  in  die 
AO  for  die  specific  window  A  as: 

p  =  Poexp(Q*|y4|) 


*  Holzmanr i  and  Cochran  (2012) 

Air  university  .  The  intellectual  and  Leadership  CenteroT  the  Air  Force 
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Literature  Review 

The  A  FIT  of  Today  is  the  Air  Force  of  Tomorrow f 


•  Extended  Spatial  Poisson  Process  Example 


r*- 


0.095  0.284  0.346 
0  0.152  0.366 
0  0  0.243 

0  0  0 
0  0  0 
0  0  0 


0.205  0.063  0.007 
0.33  0.132  0.02 

0.439  0.264  0.053 
0.39  0.469  0.141 
0  0.624  0  3  76 
0  0  1 


=  [o»5 


0184  0346  0105  0.063  0.007 


] 


*  Holzmann  and  Cochran  (2012) 

Air  University.  The  intellectual  and  Leadership  Center  or  the  Air  Force 
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Literature  Review 

1  TheAFIT  of  Today  is  the  Air  Force  of  Tomorrow* 


Matrix  Exponentiation 

•  MATLAB  &  SAS  quick...  but  not  readily  accessible. 

•  Numerically 

"  (<?  *  Ml)4 


_  V 

k-  0 


k\ 


Diagonalization 

eQ  =  S  ■  eQ  •  5_1  S 


0^-n. 


— 1 


>  Strang  (2 006 ) 


Atr  university:  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
Aim  High.  .Fly- Fight- Wm 


19 


W  Methodology  fMFIT 

The  AFIT  of  Todayis  the  Air  Force  of  Tomorrow 

•  Assumptions 

•  Methods 

•  Parameters 

•  Flow  Diagram 

•  Model 

•  Outputs 


Air  University.  The  Intellectual  and  Leadership  Center  of  the  Air  Force 
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Methodology  -  Assumptions  <-AF 1 1 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow^^^^^^^^ 

•  The  baseline  for  an  Operations  Area  (AO)  is  1  km3 

•  #  and  type  of  birds  known  to  an  acceptable  level. 

•  The  velocity  of  birds  in  the  AO  is  negligible  compared  to 
the  speed  of  the  aircraft  -  bird  is  stationary  wrt  aircraft. 

•  The  number  of  birds  in  the  AO  is  constant. 

•  An  aircraft/bird  encounter  within  the  ellipse  formed  by  the 
height  and  width  of  an  aircraft  will  be  considered  a  strike. 

•  Altitude  distributions  of  birds  (if  known)  will  following  a 
gamma  distribution  (birds  within  an  altitude  band  will  be 
distributed  via  a  Poisson  process). 

•  Small  Bird  Estimations  are  utilized  as  needed 

A)r  university.  The  intenectuai  ana  Leadership  Center  oT  the  Air  Force 
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W  Methodology- Parameters  <-AFTl 

^ ■  The  AFITof  Today  is  the  Air  Force  of 

Parameters 


2dcmi£cr 

^cce^ces 

Used  is 

L 

AO  Lcnjtfc  is  meters 

[V]  ->AO  Volume 

W 

AO  Width  s as  merer  j 

V  *^AO  Vekuae 

XI 

AO  #*  er^'r*  is  me  ten 

V  ^A0  Volume 

A_\k 

.\ircrsi  Al&aide  N5sx 

H  ^Modified  AO  Volume 

A 

.kir  cr cr  Alcaide  N&a 

%  ^  Me  diced  AO  Vefcuoe 

A_W 

Airerss 

A  Aircraft  Voumc 

aju 

Air  eras 

A  Aircraft  Vacate 

i-xreer  .Ares 

.kirerss  FrottsI  *A 

A  Aircraft  Vdurae 

A Jb 

Aker  aft  Btitesce 

A  Aucnft  Vduac 
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Tael  •  srei  js  AO 

7. 4  ktaisn  rancsen 

Bc _ c 

Gamma  Alyha  -.sltjc 

A  caMAl*  ktsssrr.  Funcccn  me  if  tea  veUsne 

a 

Cremms  Bets  vsiwe 

X’  Ate  istessir.  Fusecec  medfisd  rolomc 

B„7 

W  cf  7 

A  jseensin  Funcccn  bird  a  Iscas 

ip  ceres  3Usk  ka: 

*  toes',  risk 'r-.  jjeeeei 

Calculations 

IVI  =  L  *  W  *  H  {length  *  width  *  height} 

Ivl  =  L  *  W  *  {A_Mx  -  A_Mn)  {length  *  width  *  attitude  band  height} 

IAI  =  '4{A_W)  *  */;{A_H)*  )*tt*(A_D}  {area  of  eftpse  *  distance} 

A=T  m 

A  =  .T  *T)/M 

A  k—A--  =  [(BvaT  *T  )*(GafnmaCDF{A_Mx}-GammaCDF{A_Mn)]/M 
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Methodology  -  Model 


1  TheAFIT  of  Today  is  die  Air  Force  of  Tomorrow* 


Case#1 

1)  Calculate  hit  probab*ty 

(x-iai r 

Pn  =  e  x  IA|  * - — - 


Case  #2 

1)  Calculate  hit  probabdty  per  specific  speces 


Pspeoes  (n)  =  e  'W  * 


sp«oe  5  *  |A|)” 


n! 


X  (assented  nasa  percesEae  of  specific  bird  in  the  AD) 


At  3  cases  need  AC  dat3  &  Ajo^f^dara. inputs 
Fotowng  are  the  extra  mputs  needed  by  case 
and  the  respective  output:: 

Case  1 1nput  =  Total  beds 
Case  1  Output  =  Probab*ty  of  Hit 

Case  2  Inputs  =  Total  beds  &  Spec  es  % 

Case  2  Output  =  Hit  probabdity  and  Risk  Score 

Case 3  Inputs  =  Total  beds.  Speces  %.  Gamma 
Ostrixitkxi  and  Aircraft  altitude  band 
Case  3  Output  =  Hit  ptobabduy  and  Risk  Score 
(fideity  on  the  altitude  band) 


2)  Aggrega:e->  number  of  beds  xprobabtty  of  hit 
(efcmmat mg  double  counts) x  speces  specific  risk  level 

N 

R-species  —  ^  *  Pspedes  (*0  *  species 

n— 0 

3)  Sum  ove' al  speces 


R  ^  Rspecies 


AD  species 


Air  university:  m«  intellectual  and  leadership  Center  of  the  Air  Force 

Aim High.  .Fly- Fight- Win  25 


71 


W  Methodology- Inputs 

The  AFIT  of  Today  is  the  Air  Force  of  TomorrowF^^ 


CASE#1  -  Inputs 


Calculate 


Aw  University  The  intellectual  and  leadership  Center  of  the  Air  Force 
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W  Methodology- Outputs  <-AH  1 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow^^^"^^^^^- 

CASE#1  -  Output 
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The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow^^^^^^^^^- 

CASE#2  -  Output 


Probability  of  No  hit 

0.89 

Prob  of  at  least  1  hit 

0.11 
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Analysis 
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Risk  Score  (Vulture) 
Altitude  Band  (300  -  600  ml 
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Conclusion 

TheAFIT  of  Today  is  the  Air  Force  of  Tomorrow! 


•  The  Spatial  Poisson  Process  accounts  for  the  spatial 
distribution  of  birds  within  a  bounded  area  of  operations. 

•  The  Extended  Spatial  Poisson  Process  allows  for  over¬ 
dispersed  and  under-dispersed  object  densities  in  this 
bounded  space  (the  AO). 

•  The  Relative  Risk  Scores  account  for  differing  hazards 
associated  with  striking  different  species  of  birds. 

•  The  Gamma  Altitude  Distributions  provides  further  fidelity 
when  determining  the  hazard  associated  with  a  particular 
altitude  band  to  be  flown. 

•  The  model  allows  for  varying  levels  of  user  input  and  provides 
numerous  outputs  to  include  collision  probabilities,  a  raw  risk 
score,  a  standardized  and  graphical  relative  risk  score,  a 
gamma  distribution  determination  tool  and  a  Risk  Filtering  and 
Ranking  Method  graphic. 
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Areas  for  Future  Research  •  .-Al  l  I 

The  AFIT  of  Today  is  the  Air  Force  of  Tomorrow^^^^^^^^^ 

•  Account  forthe  varying  flock  sizes 

•  Shot  noise  Cox  processes  and  the  finite  Gibbs  point 
process  could  be  utilized  to  model  the  clustering  of  birds 
within  an  AO. 

•  Flock  shapes  and  density  distributions  within  the  flock 
could  be  studied.  We  treat  them  as  random,  but  there 
may  be  structure. 


*  l.ialier  &  Waagepa!  arson.  1 2006) 
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